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ABSTRACT 

A  substantial  amount  of  work  into  the  statistical  analysis  of  both 
visual  and  infra-red  imagery  has  been  undertaken  in  the  recent  years 
at  DSTO.  This  report  is  a  summary  of  the  results  from  the  analysis  of 
a  trial  held  in  Northern  Australia  in  1997.  The  first  and  second  order 
statistics  of  land  background  infra-red  pixel  intensities  were 
estimated  and  Gaussian  histogram  and  power-law  auto-covariance 
models  were  rigorously  tested,  before  a  stochastic  simulation  of  the 
data  was  performed.  As  the  simulation  technique  required  Gaussian 
statistics,  a  point-wise  non-linear  transform  was  applied  to  those 
ensemble  images  with  non-Gaussian  histograms. 
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Statistical  Analysis  of  Northern  Australian  Land 

Backgrounds 


Executive  Summary 


Knowledge  of  the  background  statistics  is  required  for  the  design  of  suitable 
background  clutter  suppression  and  target  extraction  filters.  This  information  can  also 
be  used  for  computer  generation  of  realistic  stochastic  models  for  system  performance 
analysis  and  generation  of  large  numbers  of  random  backgrounds  for  Monte-Carlo 
testing.  This  report  is  a  summary  of  the  first  and  second  order  statistics  of  natural 
terrain  imagery  in  Northern  Australia  during  the  1997  trial.  The  statistics  are  also 
applied  to  generate  synthetic  images  over  these  terrains. 

The  ensembles  were  formed  by  parsing  each  recorded  rim  and  grouping  images 
containing  the  predominate  vegetation  class.  Using  the  assumption  that  the  underlying 
random  processes  are  stationary  in  the  wide-sense,  the  ensemble  histograms  and  auto¬ 
covariance  can  be  calculated.  These  then  become  estimates  of  the  true  distribution  and 
correlation  for  that  vegetation  classes. 

Using  Dietrich  and  Newsam's,  [1]  technique  for  generating  arbitrary  stationary 
Gaussian  random  fields  using  the  first  and  second  order  statistics  of  an  existing  image, 
synthetic  images  can  be  simulated.  For  sequences  without  an  underlying  Gaussian 
distribution,  the  non-linear  transform  developed  by  Chappie  et  al.  [2]  can  be  used  to 
generate  a  non-Gaussian  random  field  with  a  marginal  probability  density  and 
correlation  structure  that  agrees  with  that  of  actual  IR  background  imagery. 

The  analysis  shows  that  the  majority  of  natural  terrain  background  is  non-Gaussian 
and  the  power-law  model  for  the  auto-covariance  matches  for  the  majority  of  the 
vegetation  classes.  Most  of  the  transformed  synthetic  images  appear  to  visually  match 
the  sample  images,  but  there  does  appear  to  be  some  difference  in  the  brightness  for 
images  which  are  not  transformed  but  pass  the  Gaussian  hypothesis  test.  However,  if 
the  transformation  is  applied  to  these  images  before  synthesis,  reasonable  simulated 
images  result. 
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1.  Introduction 


1.1  Aim 

Many  advanced  image  processing  algorithms  used  for  automatic  detection/ recognition 
of  objects  in  natural  terrain,  depend  on  the  spatial  statistics  of  the  background,  [3]. 
Knowledge  of  the  background  statistics  is  required  for  the  design  of  suitable 
background  clutter  suppression  and  target  extraction  filters,  [4].  This  information  can 
also  be  used  for  computer  generation  of  realistic  stochastic  models  for  system 
performance  analysis  and  generation  of  large  numbers  of  random  backgrounds  for 
Monte-Carlo  testing,  [l]-[2]. 

To  aid  in  this  analysis,  visual  and  infra-red  image  data  in  the  0.4pm  to  0.8pm  (visible), 
3pm  to  5pm  (MWIR)  and  8pm  to  12pm  (LWIR)  bands  were  collected  during  a  series  of 
trials  in  both  Southern  and  Northern  Australia.  The  visual  trial  was  conducted  in  July 
1993,  followed  by  four  trials  using  infra-red  bands  in  December  1995,  February  1996, 
July  1997  and  August  1998.  The  first  two  of  these  trials  were  conducted  by  the  Wide 
Area  Surveillance  Systems  Division  (WASSD),  now  Surveillance  Systems  Division 
(SSD)  and  the  other  three  by  the  Weapons  Systems  Division  (WSD). 

A  total  of  eight  reports,  [l]-[2]  and  [4]-[8],  have  been  written  focussing  on  different 
parts  of  this  analysis,  each  one  using  only  a  small  subset  of  the  data  which  has  been 
collected.  This  report  aims  to  fill  in  the  gaps,  by  analysing  an  entire  data  set  from  the 
1997  trial  and  rigorously  test  the  main  findings  of  papers  [l]-[2]  and  [5]-[8],  The  other 
two  reports  focus  on  specific  application  of  the  results  and  will  not  be  considered  here. 


1.2  Trials  Background 

The  visual  trial  in  1993  and  the  MWIR  trial  in  1995  were  conducted  by  the  Image 
Exploitation  group  in  WASSD.  Analyses  of  the  data  from  these  trials  led  to  papers  [1], 
[5]  and  [6].  This  work  was  then  continued  by  both  the  Electro-Optical  Seekers  and 
Missile  Simulation  groups  in  WSD.  A  further  three  trials  were  conducted  in  the  MWIR 
and  LWIR  in  1996,  1997  and  1998,  which  led  to  a  further  two  papers,  [7]-[8J.  More 
recently,  the  work  in  [1]  has  been  extended  and  papers  [2]- [4]  have  been  written. 

In  the  WSD  trials,  a  pair  of  Agema  900  imaging  radiometers  were  used  to  capture  the 
imagery  in  both  the  MWIR  and  the  LWIR  bands.  The  imagers  were  mounted  inside  a 
RAAF  Dakota  Aircraft,  pointing  sideways  out  of  the  parachute  door.  The  imagers'  lines 
of  sights  (LOS)  were  aligned,  but  not  'pixel  registered'1,  due  to  different  aberrations  in 
their  optical  systems.  During  data  collection,  the  aircraft  travelled  at  a  speed  of 
approximately  120  knots. 


1. 


If  the  imagers  were  'pixel  registered',  then  every  pixel  recorded  by  the  LWIR  system 
has  exactly  one  corresponding  pixel  in  the  same  row  and  column  recorded  by  the 
MWIR  and  vice  versa. 
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The  purpose  of  the  initial  WSD  trial  in  February  1996  was  to  collect  data  at  various 
sensor  look-down  angles  between  10°  and  60°.  It  turned  out  that  the  range  of 
depression  angles  between  25°  and  50°  could  not  be  used,  as  the  imagers'  field  of  views 
(FOV)  were  affected  by  one  of  the  aircraft  engine's  exhaust  gasses.  Second,  data 
collection  at  shallow  look  down  angles  smaller  then  15°  resulted  in  unusable  images 
due  to  high  atmospheric  attenuation.  Data  collection  during  the  1996  wet  season  trial 
was  therefore  limited  to  sensor  look-down  angles  of  60°. 

As  the  images  recorded  in  the  1996  trials  were  generally  of  low  contrast  and  very  noisy, 
due  to  high  atmospheric  absorption  and  low  temperature  contrast  due  to  cloud  cover, 
the  trial  was  repeated  in  July  1997  to  collect  dry  season  background  data  of  the  same 
areas  overflown  in  February  1996.  The  image  data  collected  in  the  dry  season  proved  to 
be  of  far  better  quality  due  to  lower  atmospheric  absorption  and  higher  temperature 
contrast  caused  by  high  solar  loading.  The  third  trial  in  August  1998  collected  data  at 
shallow  sensor  depression  angles  around  15°.  In  addition,  any  incomplete  data  from 
previous  data  sets  were  collected.  The  exact  locations  of  the  trials  are  given  in  the 
classified  appendix. 

The  geographical  area  over  which  image  data  was  collected,  was  selected  from  a 
vegetation  map  of  the  Northern  Territory  (NT)  [see  classified  appendix].  This  approach 
was  used  to  ensure  that  the  vegetation  types  most  commonly  found  in  the  NT  were 
covered. 

1.3  Imager  Details 

The  Agema  900  uses  a  pair  of  scanned  cooled  detectors  made  of  Indium-Antimonide 
(InSb)  for  the  MWIR  and  Cadmium-Mercury-Telluride  (CMT)  for  the  LWIR.  Image 
frames  are  collected  at  a  rate  of  2Hz  with  a  FOV  of  2.5°  vertical  by  5°  horizontal. 

As  the  system's  modulation  transfer  function  affects  the  horizontal  spatial  resolution, 
each  frame  was  recorded  digitally  as  a  136x272  pixel  image  to  maintain  the  correct 
aspect  ratio.  This  gave  a  collection  rate  with  approximately  50%  overlap  between  one 
image  and  the  next  and  modified  the  frame  rate  to  approximately  1Hz. 

The  radiance  distribution  projected  onto  the  imaging  plane  of  the  sensor  is  related  to 
the  radiance  distribution  from  the  terrain  (self-emitted  plus  reflected  solar)  by 
convolution  with  the  optics  point-spread  function.  For  the  purpose  of  the  analysis 
contained  within  this  report,  the  effect  of  this  convolution  is  insignificant,  and  has  been 
neglected. 

Each  pixel  has  a  12-bit  dynamic  range  before  the  AGEMA  software  linearises  the 
response  on  the  detector.  A  non-linear  mapping  from  the  integer  range  0-4095  (12  bit) 
to  0-8191  (13-bit)  causes  a  large  number  of  empty  bins  in  the  ensemble  histogram. 
These  empty  bins  are  later  removed  by  a  median  filter  of  length  five.  Although  not 
mathematically  exact,  the  effect  on  the  intensity  distribution  is  negligible.  After 
linearisation,  the  measured  pixel  intensities  are  taken  to  be  directly  related  to  the 
intensity  of  the  radiation  emanating  from  the  scene,  including  atmospheric 
propagation  effects.  No  attempt  was  made  to  account  for  atmosphetic  attenuation  and 
path  radiance  as  it  does  not  significantly  affect  the  shape  of  either  the  histogram  or  the 
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auto-covariance  function. 


1.4  Document  Overview 

Most  current  generation  image  processing  algorithms  implemented  in  electro-optic 
surveillance  systems  use  only  first  and  second  order  statistics.  Therefore,  the  analysis 
of  the  data  in  this  report  has  been  limited  to  the  first  and  second  order  statistics  of  the 
radiance  intensity  distribution. 

Previous  results  indicate  that  the  radiance  intensity  distributions  of  the  natural  terrain 
backgrounds  from  the  WSD  trials1  are  nearly  always  wide-sense  stationary  but  not 
necessarily  Gaussian,  [7]-[8].  Therefore,  the  analysis  in  this  report  starts  by  assuming 
the  ensemble  images  are  wide-sense  stationary.  The  ensemble  histogram  is  then 
formed  and  the  hypothesis  that  the  distribution  is  Gaussian  is  tested.  The  hypothesis 
threshold  is  tested  by  the  Root  Mean  Square  (RMS)  error  and  a  Chi-squared  test. 

The  auto-covariance  estimate  is  split  into  horizontal  and  vertical  slices.  A  power-law 
model  is  applied  to  the  slices  and  a  threshold  test  based  on  the  Root  Mean  Square 
(RMS)  error  is  used  to  test  the  model.  A  full  auto-covariance  calculation  of  a  single 
sample  image  is  then  used  to  create  a  stochastic  simulation.  For  images  with  non- 
Gaussian  statistics,  a  point-wise  non-linear  transformation  is  applied. 

The  first  part  of  this  report  explains  the  background  to  this  work  and  a  detailed 
description  of  the  different  analysis  techniques.  The  second  part  is  a  listing  of  the 
results  and  comparisons  between  MWIR  and  LWIR  histograms,  auto-covariance  slices 
and  stochastic  simulations  for  each  different  vegetation  class. 


2.  Analysis  Techniques 

2.1  Forming  an  Ensemble 

Given  the  imagers'  FOVs  (2.5°  vertical  x  5°  horizontal),  the  number  of  pixels  (136 
vertically  x  272  horizontally)  and  a  known  look-down  angle  and  height,  the  area 
covered  by  a  single  pixel  on  the  ground  can  be  calculated  for  each  run.  The  pixel  size 
strictly  applies  to  pixels  closest  to  the  image  centre  only.  For  pixels  on  the  top  and 
bottom  rows,  there  is  a  variation  of  approximately  +4%  that  was  not  corrected  for. 

In  this  analysis,  each  recorded  run  was  parsed  and  images  of  the  predominate 
vegetation  class  were  extracted  to  form  an  ensemble.  Runs  with  the  same  predominate 
vegetation  class  and  pixel  size,  were  combined  to  form  a  larger  ensemble.  This 
approach  lets  each  vegetation  class  be  independent. 


1.  Wide-sense  stationarity  depends  on  many  things,  such  as  sensor  orientation  and  the  field 
of  view.  These  conditions  were  met  for  the  WSD  trial. 
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2.2  Random  Field  Analysis 

In  the  following,  each  image  in  an  ensemble  is  treated  as  a  truncated  segment  of  a 
different  realisation  of  a  two-dimensional  discrete-space  random  field.  With  this 
assumption, 

x‘(m,n)  (1) 

with  m  =  1,2,... ,136;  n  =  1,2,... ,272  and  i  =  1,2,... ,L 

where  L  is  the  number  of  images  in  the  ensemble,  represents  the  pixel  at  (m,n)  for  the  i- 
th  image. 

This  approach  is  different  from  the  approach  used  by  Yaglom  [9],  who  assumes  in  his 
analysis  that  every  image  in  the  ensemble  is  a  single  realisation  of  an  ergodic  random 
field.  Although,  the  ergodic  assumption  is  intuitively  appealing  as  the  sensors  record 
different  patches  of  die  same  class  of  terrain,  difficulties  will  arise  later  if  the 
assumption  does  not  hold.  The  approach  used  for  this  analysis  does  not  rely  on  this 
assumption,  therefore  the  results  will  remain  valid  regardless  of  whether  the 
underlying  process  is  ergodic. 

Two  different  approaches  can  also  be  used  to  undertake  the  correlation  analysis  on  the 
image  ensemble.  First,  correlation  analysis  can  be  carried  out  after  subtracting  from 
every  individual  image  in  the  ensemble  its  own  mean,  [10],  and  second  the  ensemble 
mean  can  be  subtracted  from  the  individual  images,  [2]  and  [6].  The  first  approach 
makes  the  inherent  assumption  that  the  auto-covariance  estimates  decay  to 
insignificantly  small  values  within  the  spatial  extent  of  one  image  frame.  As  this 
requirement  cannot  be  guaranteed  for  the  type  of  imagery  analysed  here,  the  second 
approach  is  being  taken  in  this  study.  If  the  auto-covariance  estimate  indeed  drops  to 
an  insignificant  level  within  the  spatial  extent  of  one  image  frame,  then  both 
approaches  will  give  identical  results.  If,  on  the  other  hand,  the  correlation  structures 
in  the  image  extend  beyond  one  image  frame,  the  approach  used  here  will  give  a  better 
approximation  of  the  auto-covariance  estimate.  As  the  spatial  extent  of  the  ensembles' 
correlation  was  not  known  at  the  beginning  of  this  study,  the  second  approach  is  used. 


2.3  The  Gaussian  Hypothesis 

2.3.1  Background 

The  analysis  of  the  first  and  second  order  statistics  of  the  image  ensembles  is  based  on 
the  assumption  that  the  underlying  random  processes  are  stationary  in  the  wide-sense. 
Mathematically,  the  mean  calculated  over  an  entire  ensemble  of  infinitely  many 
realisations  of  the  random  field, 

(2) 

/u  =  E[x(m,n)]r 

and  its  auto-covariance  at  lag  ( k,l ) , 
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C(k,  l )  =  E[(x(m,  n )  -  / i)(x(m  +  k, «  +  /)-//)] ,  O' 

are  independent  of  position  (m,n). 

In  this  report,  the  image  ensembles  were  not  tested  for  wide-sense  stationarity  due  to 
the  findings  in  [7]-[8].  In  addition,  if  the  random  processes  were  Gaussian  it  would 
imply  strict-sense  stationarity.  This  condition  is  tested  rigorously  in  this  report  by 
using  an  RMS  comparison  and  a  Chi-squared  test. 


2.3.2  Histogram 

To  test  the  Gaussian  hypothesis,  an  ensemble  histogram  is  calculated  and  overlayed 
with  a  Gaussian  plot  as  shown  in  Figure  1.  The  ensemble  histogram  is  obtained  by 
grouping  all  of  the  pixels  from  each  of  the  images  in  the  ensemble,  and  then  counting 
the  fraction  of  the  pixels  that  have  grey-level  0,1,2,. .  .,8191. 

The  Gaussian  model  is  based  on  a  single-variate  distribution  with  the  mean  (//)  and 

variance  (a2 ),  which  are  determined  by  minimising  the  RMS  between  the  Gaussian 
model  and  the  ensemble  histogram.  The  Gaussian  model  is  given  by 


h2  (x)  = 


_1 _ 

yjlncr2 


exp 


(*~/Q2 

2  a2 


(4) 


The  RMS  error  can  then  be  calculated  by  the  following: 


I  1 

where  hj(x)  is  the  ensemble  histogram  and  hi(x)  is  the  Gaussian  model  at  position  x. 
Since  both  hi(x)  and  hi(x)  are  normalised,  the  RMS  value  is  expressed  as  a  percentage. 

The  Chi-squared  test  needs  to  be  constructed  with  bins  of  at  least  five  samples 
(Cochran's  Rule,  [11]).  The  ensemble  histogram  is  hence  trimmed  to  remove  any  points 

smaller  then  five  and  the  result  is  a  new  histogram  /*,  with  its  bins  now  labelled  m(m  = 
1,..  ,,M).  The  trimmed  points  are  not  used  in  the  comparison. 

The  test  statistic  is  then  computed  as: 

.  _  y  Pi  P0  ~  h2  (m))  ^ 

^  h2(m) 

with  corresponding  P-value:  P(zL- 1-2  >  T)  ■  The  degrees  of  freedom  (DOF)  is  M-l-2,  as 
two  parameters  were  estimated  in  the  Gaussian  model. 
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Histogram 


Figure  1  -  Histogram  with  Gaussian  overlay 
2.3.3  Determining  Threshold 

The  hypothesis  threshold  for  whether  the  ensemble  histogram  is  Gaussian  or  not  is 
determined  by  the  RMS  value  and  the  Chi-squared  value.  The  limits  imposed  were  less 
than  10%  error  for  the  RMS  and  a  Chi-squared  value  greater  than  5%. 

The  results  of  this  test  were  also  used  to  determine  if  the  image  needed  to  be 
transformed  when  constructing  a  stochastic  simulation. 


2.4  First  and  Second  Order  Statistics 

The  mean  and  auto-covariance  of  a  wide-sense  stationary  field  can  theoretically  be 
obtained  using  (2)  and  (3).  However,  since  they  require  an  infinite  number  of 
realisations  of  the  random  processes,  these  equations  must  be  modified  to  reflect  the 
fact  that  the  sample  size  is  finite. 

Given  that  the  underlying  random  processes  are  wide-sense  stationary  and  that  the 
individual  images  are  statistically  independent,  the  best  possible  estimate  of  the  mean 
of  the  random  process  is  given  by 


1 


L  M  N 


x'(m,n), 


(7) 
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with  variance 

<r?=£[oi-/,)j]«I^h_|;ai'-1a)2,  (8) 

where  jj!  is  the  mean  of  the  i-th  image. 

This  estimate  of  the  mean  approaches  the  true  mean,  fi,  of  the  random  field  as  the  size 
of  the  ensemble  approaches  infinity. 

An  estimate  of  the  auto-covariance  function  is  obtained  from  the  finite  sample  size  in  a 
similar  way 


C(k,l) 


1 


L  M-kN-l 


S  S Z (x‘ (m’ ") _ £)(*' (*n  +  k,n  +  l)~  ft), 


(9) 


with  variance 


<7?  =  Ei(c(k,i)-c(k,i)Y]  *  1  -il(c‘(k,r)-c(k,D)2 , 

L*(L- 1)^ 


(10) 


where  C‘  ( k ,  l )  is  the  auto-covariance  of  the  z-th  image. 

Again  as  the  sample  size  increases  to  infinity,  this  estimate  approaches  the  true  auto¬ 
covariance  function  C(k,l). 

As  the  estimates  C(k,l )  and  //  are  obtained  from  a  large  amount  of  independent 
images,  it  follows  from  the  Central  Limit  Theorem,  [12]  that  both  can  be  taken  to  be 
Gaussian  random  variables  with  mean  //  and  auto-covariance  C(k,l),  respectively. 
Figure  2  is  an  example  of  the  auto-correlation  function,  plotted  in  three  dimensions. 
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3D  Auto-covariance  of  Image  Ensemble 


As  the  lag  increases  away  from  the  central  peak,  there  are  increasing  fluctuations. 
These  fluctuations  represent  increasing  statistical  variation  and  is  associated  with  the 
fact  that  fewer  data  samples  are  available  at  larger  lags.  To  simplify  the  analysis,  the 
auto-covariance  is  calculated  for  each  horizontal  and  vertical  lag.  They  are  then 
averaged  to  give  horizontal  and  vertical  slices.  Figure  3  demonstrates  this  process  with 
one  standard  deviation  error  bars  also  shown. 


Auto-covariance  -  Horizontal  Auto-covariance  -  Vertical 


Figure  3  -  Auto-covariance  for  Horizontal  and  Vertical  slices 
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2.4.1  Auto-covariance  Models 

In  previously  papers,  [5]-[7],[10]  attempts  were  made  to  model  the  auto-covariance 
functions  both  as  an  exponential  and  with  a  power-law  model. 

The  exponential  function  was  suggested  by  N.  Ben-Yosef  et  al.  [10]  and  has  an  equation 
of  the  form 


C(i.0  =  ^*exp(-^(A,t)2  +(A,/)2),  <n> 

where  Ax  and  Ay  are  the  horizontal  and  vertical  dimensions  of  a  pixel  and  A  >  0  and  y 
>  0  are  constants. 


It  was  found  in  [5]-[6]  that  this  model  did  not  match  the  auto-covariance  estimates 
measured.  Instead,  a  power-law  model  was  tested  by  [5]  and  [7]  for  an  isotropic 
random  field,  which  exhibits  a  power-law  fall-off 


C(k,l) 


a 

C b2  +  k2+i2y ' 


(12) 


where  a  >  0,b>  0  and  v>  0  are  constants. 

It  was  found  in  [5]  and  [7],  that  the  power-law  rather  then  the  exponential  model 
provided  a  better  match  for  the  auto-covariance  estimates  measured.  Consequently, 
this  is  the  model  that  we  have  decided  to  use.  The  method  for  matching  the  power-law 
model  to  the  data  was  simplified  by  reducing  the  number  of  parameters  by  one.  The 
horizontal  and  vertical  slices  were  used  to  test  the  simplified  model,  which  now 
reduces  to 


C(k)  = 


M(b2  +k2)v 


and  C(l)  = 


2\v 


N(b  +n 


(13) 


The  best  fit  parameters  were  obtained  by  a  numeric  maximum  likelihood  estimation 
technique  on  the  correlation  data  over  about  half  of  the  available  range  of  lags,  [13]. 
They  were  then  adjusted  to  minimize  the  RMS  between  the  auto-covariance  estimate 
and  the  model.  Figure  4  is  an  example  of  two  power-law  fits. 

2.4.2  Determining  Threshold 

A  basic  threshold  test  was  used  for  determing  if  the  power  law  model  fit  the  data.  The 
RMS  value  is  3%  and  was  chosen  by  trial  and  error. 
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Auto-covariance  -  Horizontal  Auto-covariance  -  Vertical 


Figure  4  -  Power-law  model  for  Horizontal  and  Vertical  Auto-covariance  slices 

2.5  Stochastic  Simulation 


2.5.1  Background 

Stochastic  image  models  play  an  important  role  in  many  areas  of  statistical  image 
processing,  [14].  Dietrich  and  Newsam,  [1]  investigated  this  and  developed  a  technique 
for  generating  arbitrary  stationary  Gaussian  random  fields  using  the  first  and  second 
order  statistics  of  an  existing  image.  This  technique  was  used  to  generate  synthetic 
images  of  the  visual  data  collected  from  the  1993  trial  and  later  in  the  1995  MWIR  trial, 
[6]. 

More  recently,  a  non-linear  transform  was  developed  by  Chappie  et  al.  [2]  to  generate  a 
non-Gaussian  random  field  with  a  marginal  probability  density  and  correlation 
structure  that  agrees  with  that  of  actual  IR  background  imagery.  This  technique  was 
then  used  for  target  detection/ recognition  in  [3]-[4]. 


2.5.2  Stochastic  Image  Simulation  of  Gaussian  Random  Fields 

The  simulation  process  begins  by  calculating  the  auto-covariance  estimate  for  an 
ensemble  image.  It  is  then  multiplied  by  a  circular  super-Gaussian  taper  to  reduce 
artefacts  associated  with  the  variance  in  the  estimate  at  large  lags.  This  is  then  used  as 
the  correlation  structure  in  the  next  part  of  the  simulation.  The  circular  super-Gaussian 
taper,  which  acts  as  a  window  to  the  auto-covariance  estimate,  is  calculated  by 


w(k )  =  exp 


1*1 


(14) 
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where  k  is  the  auto-correlation  and  P>  0  represents  the  radius  of  the  taper. 

In  this  analysis,  a  value  of  P  was  chosen  to  include  most  of  the  slowly  falling 
correlation  tail.  The  process  then  continues  by  computing  the  Fast  Fourier  Transform 
(FFT)  of  the  auto-covariance  estimate,  setting  any  negative  values  to  zero,  performing 
an  element-wise  square  root,  and  then  multiplying  the  result  by  a  white  noise  random 
function  with  Gaussian  statistics.  The  FFT  of  the  resulting  matrix  is  then  separated  into 
real  and  imaginary  parts,  both  of  which  are  independent  realisations  of  the  stationary 
zero-mean  Gaussian  random  field  with  the  required  correlation  structure.  The 
simulation  is  completed  by  adding  the  estimate  of  the  mean  of  the  random  field.  Figure 
5  is  a  visual  representation  of  this  process. 


Irnag.  Part  +  Mean 


Figure  5  -  Gaussian  Simulation  Process 

There  are  two  sources  of  error  in  this  simulation  approach.  The  first  is  due  to  the 
tapering  of  the  correlation  data  and  is  unavoidable  because  of  the  fact  that  the  variance 
in  an  auto-covariance  estimate  will  always  increase  at  large  lags  as  a  result  of  the  small 
number  of  data  samples.  The  second  is  due  to  setting  to  zero  the  negative  values  in  the 
FFT  of  the  tapered  auto-covariance  estimate.  Bertilone  et  al.  [6]  has  shown  the  latter 
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error  has  a  magnitude  of  lO6,  which  is  insignificant  in  this  simulation.  Figure  6  is  a 
comparison  between  a  sample  image  and  a  simulated  one. 


Sample  image  Simulated  image 


50  100  150  200  250  50  100  150  200  250 


Figure  6  -  Comparison  of  Sample  and  Simulated  Images 
2.5.3  Non-Linear  Transform 


One  of  the  conditions  for  Dietrich  and  Newsam's  method  to  work  is  for  the  ensemble 
image  to  have  Gaussian  statistics.  If  the  simulation  is  performed  on  a  sample  image 
without  this  distribution,  many  features  are  visually  incorrect  and  the  histogram  of  the 
new  image  is  very  different  from  the  sample  image.  Chappie  et  al.  [2]  has  developed  a 
technique  whereby  the  images  are  transformed  by  a  point-wise  non-linear  transform  to 
generate  an  ensemble  with  a  Gaussian  histogram.  After  the  image  has  been  simulated, 
an  inverse  non-linear  transform  is  used  to  restore  the  original  histogram.  Figure  7 
shows  a  comparison  between  a  synthesized  image  with  and  without  transformation.  It 
is  clear  in  Figure  7b  that  without  the  transform,  the  contrast  is  lacking  and  much  of  the 
detail  is  not  represented  correctly. 


Synthesized  Image  (Transformed)  Synthesized  Image  (Not  Transformed) 

20 
40 

60 

80 
100 
120 

50  100  150  200  250 


(a) 


(b) 


Figure  7  -  Comparison  of  Simulated  Image  (Transformed  vs.  Not  Transformed) 
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The  transformation  function  F(x)  is  used  to  convert  values  X(m,n)  to  values  Y(m,n)  = 
F(X(m,n))  with  a  Gaussian  distribution  with  zero  mean  and  a  variance  of  one.  For  a 
probability  density  Px  (x),  the  required  transformation  is 


F(x)  =  Serf-'  (2  j ‘Px  (x')dx'  - 1) .  (15) 

o 

A  new  image  is  then  created  with  this  distribution  and  is  substituted  for  the  ensemble 
image.  To  recreate  the  final  images,  the  inverse  non-linear  transform,  F~ 1  (x)  is  applied 
to  the  simulated  images.  Figure  8  shows  typical  plots  of  the  non-linear  transform  and 
the  inverse  non-linear  transform. 


Non-linear  Transform  Non-linear  Inverse-Transform 


Figure  8  -  Transform  and  Inverse  Transform 


For  a  sample  ensemble  image  with  a  non-Gaussian  distribution.  Figure  9  shows  the 
histogram  after  using  the  above  transform  to  create  a  sample  image. 
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Histogram  of  sample  image  X10'3  Histogram  of  transformed  sample  image 


Figure  9  -  Original  and  Transformed  Histograms  of  a  Sample  Image 


After  the  transformation.  Figure  9b  now  shows  a  Gaussian  distribution  with  zero  mean 
and  a  variance  of  one.  After  creating  the  simulated  images.  Figure  10  shows  a 
comparison  between  the  histogram  of  the  original  simulated  image  and  one  after  the 
inverse  transform  has  been  applied. 


x  i  o3  Histogram  of  transformed  simulated  image 


(a) 


(b) 


Figure  10-  Original  and  Transformed  Histograms  of  a  Simulated  Image 


Its  clear  that  the  original  distribution  has  been  restored  with  very  little  difference 
between  Figure  9a  and  Figure  10b.  To  verify  the  second  order  statistics.  Figure  11  is  a 
comparison  between  the  sample  and  simulated  horizontal  and  vertical  auto-covariance 
slices. 
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x-104  Auto-covariance  -  Horizontal  (Comparison)  X104  Auto-covariance  -  Vertical  (Comparison) 


Figure  11  -  Auto-covariance  Slices  of  the  Sample  and  Simulated  Image 


The  peak  at  the  start  is  well  matched  in  both  plots,  but  the  remainder  of  the  tail  is 
subject  to  a  statistical  fluctuation  associated  with  the  small  sample  size  and  windowing 
the  centre  of  the  auto-covariance  estimate. 


3.  1997  Trial 

3.1  Sequence  Summary  from  AGEMA  900  IMAGER 

During  the  1997  trial,  three  groups  of  runs  were  taken  on  different  days  and  at 
different  times  in  the  day.  The  predominate  vegetation  class  has  been  determined  and 
forms  the  basis  of  the  analysis  for  each  run  (see  section  2.1  for  information  on  how  the 
classes  were  chosen).  Table  1  is  a  brief  summary  of  each  run.  The  key  for  the  vegetation 
classes  is  in  the  classified  appendix. 
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Table  1  -  1997  Sequence  Summary 


Flight  Date  and 
Conditions 

Run 

Number 

Altitude 

(feet) 

Depression 

Angle 

(degrees) 

Vegetation 

Classes  in  Run 

Predominate 
Vegetation  Class 

FLIGHT  DATE: 
30-6-97,  NIGHT 

1 

10K 

60 

9, 105, 11, 4, 48 

4 

2 

10K 

60 

48, 4, 105 

4 

3 

5K 

60 

9, 105, 11, 4, 

9 

4 

5K 

60 

54, 48, 4 

54 

5 

1.7K 

60 

11,4 

11 

6 

1.7K 

15 

48 

48 

7 

1.0K 

tracking 

N/A 

N/A 

8 

1.0K 

tracking 

N/A 

N/A 

FLIGHT  DATE: 
1-7-97,  DAY 

9 

10K 

60 

105, 21, 11, 32, 4, 9 

32 

FLIGHT  DATE: 
3-7-97,  DAY 

10 

10K 

60 

48, 11, 21, 105 

11 

11 

10K 

60 

21, 15,  51 

21 

12 

10K 

60 

13,  51, 15, 4,  21 

15 

13 

5K 

60 

9, 21, 11, 32 

32 

14 

5K 

60 

48, 11,  51 

11 

15 

5K 

60 

21, 15,  51 

21 

16 

5K 

60 

13,  51, 15, 4,  21 

15 

17 

10K 

60 

4, 106 

106 

18 

10K 

60 

53, 4, 106 

53 

19 

5K 

60 

4,106 

106 

20 

5K 

60 

53, 4, 106 

53 

21 

10K 

60 

18, 106 

106 

22 

10K 

60 

106, 3, 47 

3 

23 

10K 

60 

3,105 

3 

24 

5K 

60 

18, 106 

106 

25 

5K 

60 

106, 3, 47 

3 

26 

5K 

60 

3,105 

3 

3.2  Explanation  of  Run  Analysis 

The  layout  of  the  following  section  is  in  four  parts: 

•  The  first  is  a  general  summary  of  the  run  and  the  ensemble. 

•  The  second  part  is  an  analysis  of  the  histogram  for  the  LWIR  and  the  MWIR. 

•  The  third  part  is  a  comparison  between  the  horizontal  and  vertical  auto¬ 
covariance  slices  for  the  LWIR  and  MWIR. 

•  The  final  part  is  a  comparison  between  sample  and  simulated  images  of  the 
ensemble. 

If  the  Gaussian  threshold  test  is  passed,  both  the  transformed  and  non-transformed 
simulated  images  are  compared.  In  addition  to  the  threshold  tests,  a  visual  test  is 


16 


DSTO-TR-1456 


added  to  each  of  the  comparisons.  The  rating  of  this  test  is  either  'yes',  'partial'  or  'no' 
depending  on  the  plots. 

3.3  Runs  1-2 

3.3.1  General  Information 

Flight  Date:  30-6-97  NIGHT 
Altitude:  lOKft 
Depression  Angle:  60  deg. 

Vegetation  Classes  in  Run:  9, 105, 11, 4, 48 
Predominate  Vegetation  Class:  4 

Ensemble  Size  is:  79 
Number  of  Valid  Images  is:  79 
Vertical  Resolution  is:  1.30m 
Horizontal  Resolution  is:  1.13m 

3.3.2  Histogram 


Histogram  for  ’LwrunOI-02' 


Histogram  for  'MwrunOI  -02’ 


0.018 

0.016 

0.014 

I"  0.012 

<D 

I  0.01 
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I  0.008 
“  0.006 
0.004 
0.002 
0 


Greylevel  Greylevel 


Figure  12  -  Gaussian  Fit  -  Histogram  of  Runs  01-02 


Histogram  Information  for  'LwrunOI-02' 

Normal  Model  Parameters: 

Mean  is:  1612.89 
Variance  is:  587.43 
RMS  error  is:  5.32% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Yes 
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Histogram  Information  for  'MwrunOl-02' 

Normal  Model  Parameters: 

Mean  is:  555.06 
Variance  is:  107.00 
RMS  error  is:  10.17% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Yes 

3.3.3  Auto-covariance  -  Horizontal 


Auto-covariance  -  Horizontal  Slice  for'LwrunOI-02' 


Auto-covariance  -  Horizontal  Slice  for  'MwrunOl-02’ 


Figure  13  -  Power  Laxv  Fit  -  Horizontal  Slice  of  Runs  01-02 


Auto-covariance  -  Horizontal  for  'LwrunOl-02' 


Mean  Standard  deviation  between  images  is:  47.91 
Power  Law  Model  Parameters: 
a_est  is:  535.04 
b_est  is:  0.00 
w_est  is:  0.04 
RMS  error  is:  1.30% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Yes 

Auto-covariance  -  Horizontal  for  'MwrunOl-02' 


Mean  Standard  deviation  between  images  is:  8.37 
Power  Law  Model  Parameters: 
a_est  is:  89.98 
b_est  is:  0.00 
w_est  is:  0.04 
RMS  error  is:  0.91% 
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Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Yes 
3.3.4  Auto-covariance  -  Vertical 


Autocovariance  -  Vertical  Slice  for  twrunOI  -02' 


Auto-covariance  -  Vertical  Slice  for  'MwrunOI-02' 


Figure  14  -  Power  Law  Fit  -  Vertical  Slice  of  Runs  01-02 

Auto-covariance  -  Vertical  for  'LwrunOl-02' 

Mean  Standard  deviation  between  images  is:  47.72 
Power  Law  Model  Parameters: 
a_est  is:  575.59 
b_est  is:  2.41 
w_est  is:  0.05 
RMS  error  is:  0.69% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Yes 


Auto-covariance  -  Vertical  for  'MwrunOI-02' 

Mean  Standard  deviation  between  images  is:  8.29 
Power  Law  Model  Parameters: 
a_est  is:  94.15 
b_est  is:  0.00 
w_est  is:  0.04 
RMS  error  is:  1.44% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Yes 
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3.3.5  Image  Synthesis: 
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Figure  15  -  Original  Images  from  Runs  01-02 

Synthesized  Image  (Transformed)  for 'LwrunOI -02' 
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Synthesized  Image  (Not  Transformed)  for  'LwrunOI-02' 
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Figure  16  -  Synthetic  Images  for  Run  01-02 


Sample  Ensemble  Image  for  'LwrunOI -02 


50  100  150  200  250 
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3.4  Run  3 

3.4.1  General  Information 

Flight  Date:  30-6-97  NIGHT 
Altitude:  5Kft 
Depression  Angle:  60  deg. 

Vegetation  Classes  in  Run:  9, 105, 11, 4 
Predominate  Vegetation  Class:  9 

Ensemble  Size  is:  250 
Number  of  Valid  Images  is:  68 
Vertical  Resolution  is:  0.65m 
Horizontal  Resolution  is:  0.56m 


3.4.2  Histogram 


Histogram  for  'Lwrun03'  Histogram  for  'Mwrun03’ 


Histogram  Information  for  'Lwrun03' 

Normal  Model  Parameters: 

Mean  is:  1667.42 
Variance  is:  585.37 
RMS  error  is:  4.78% 

Chi2  P- value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Yes 
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Histogram  Information  for  'Mwrun03' 

Normal  Model  Parameters: 

Mean  is:  621.15 
Variance  is:  109.00 
RMS  error  is:  5.10% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Yes 

3.4.3  Auto-covariance  -  Horizontal 

Autocovariance  -  Horizontal  Slice  for  'Lwrun03'  Auto-covariance  -  Horizontal  Slice  for  'Mwnjn03' 


Figure  18  -  Poxver  Laxo  Fit  -  Horizontal  Slice  of  Run  3 

Auto-covariance  -  Horizontal  for  'Lwrun03' 

Mean  Standard  deviation  between  images  is:  59.79 
Power  Law  Model  Parameters: 
a_est  is:  654.11 
b_est  is:  1.48 
w_est  is:  0.06 
RMS  error  is:  9.25% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 

Auto-covariance  -  Horizontal  for  'Mwrun03' 

Mean  Standard  deviation  between  images  is:  9.80 
Power  Law  Model  Parameters: 
a_est  is:  112.67 
b_est  is:  0.80 
w_est  is:  0.07 
RMS  error  is:  7.45% 

Model  Accepted  by  Threshold:  No 
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Model  Accepted  Visually:  No 

3.4.4  Auto-covariance  -  Vertical 


Auto-covariance  -  Vertical  Slice  for  Iwnjn03'  Auto-covariance  -  Vertical  Slice  for  'Mwrun03' 


Figure  19-  Power  Laxo  Fit  -  Vertical  Slice  of  Run  3 


Auto-covariance  -  Vertical  for  'Lwrun03' 

Mean  Standard  deviation  between  images  is:  62.34 
Power  Law  Model  Parameters: 
a_est  is:  672.35 
b_est  is:  2.29 
w_est  is:  0.06 
RMS  error  is:  1.34% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Partially 

Auto-covariance  -  Vertical  for  'Mwrun03' 

Mean  Standard  deviation  between  images  is:  9.83 
Power  Law  Model  Parameters: 
a_est  is:  137.40 
b_est  is:  8.28 
w_est  is:  0.10 
RMS  error  is:  2.95% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Partially 
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3.4.5  Image  Synthesis: 
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Figure  20  -  Original  Images  from  Run  3 

Synthesized  Image  (Transformed)  for  ’Lwrun03' 
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Synthesized  Image  (Not  Transformed)  for  'Lwrun03' 
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Figure  21  -  Synthetic  Images  for  Run  3 
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Sample  Ensemble  Image  for  'Lwrun03 
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Sample  Ensemble  Image  for  ’Mwrun03' 
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Synthesized  Image  (Transformed)  for  'Mwrun03' 


50  100  150  200  250 

Synthesized  Image  (Not  Transformed)  for  'Mwrun03' 
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3.5  Run  4 

3.5.1  General  Information 

Flight  Date:  30-6-97  NIGHT 
Altitude:  5K  ft 
Depression  Angle:  60  deg. 
Vegetation  Classes  in  Run:  48, 4, 105 
Predominate  Vegetation  Class:  54 

Ensemble  Size  is:  251 
Number  of  Valid  Images  is:  100 
Vertical  Resolution  is:  0.65m 
Horizontal  Resolution  is:  0.56m 


3.5.2  Histogram 


Histogram  for  'Lwrun04'  Histogram  for  'Mwrun04' 


Histogram  Information  for  'Lwrun04' 

Normal  Model  Parameters: 

Mean  is:  1620.83 
Variance  is:  624.96 
RMS  error  is:  7.41% 

Chi2  P- value  is:  99.91% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Partially 
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Histogram  Information  for  'Mwrun04' 

Normal  Model  Parameters: 

Mean  is:  607.98 
Variance  is:  137.03 
RMS  error  is:  6.82% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Yes 

3.5.3  Auto-covariance  -  Horizontal 

Auto-covariance  -  Horizontal  Slice  for  'LwrunCM'  Auto-covariance  -  Horizontal  Slice  for  'Mwrun04' 


Figure  23  -  Poxver  Law  Fit  -  Horizontal  Slice  of  Run  4 


Auto-covariance  -  Horizontal  for  'Lwrun04' 

Mean  Standard  deviation  between  images  is:  63.36 
Power  Law  Model  Parameters: 
a_est  is:  622.52 
b_est  is:  0.21 
w_est  is:  0.04 
RMS  error  is:  3.33% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 

Auto-covariance  -  Horizontal  for  'Mwrun04' 

Mean  Standard  deviation  between  images  is:  11.56 
Power  Law  Model  Parameters: 
a_est  is:  236.97 
b_est  is:  30.12 
w_est  is:  0.12 
RMS  error  is:  5.09% 

Model  Accepted  by  Threshold:  No 
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Model  Accepted  Visually:  Partially 

3.5.4  Auto-covariance  -  Vertical 


Auto-covariance  -  Vertical  Slice  for  'LwrunM'  Auto-covariance  -  Vertical  Slice  for  'Mwrun04' 


Figure  24  -  Power  Lazo  Fit  -  Vertical  Slice  of  Run  4 


Auto-covariance  -  Vertical  for  'Lwrun04' 

Mean  Standard  deviation  between  images  is:  62.57 
Power  Law  Model  Parameters: 
a_est  is:  613.29 
b_est  is:  0.00 
w_est  is:  0.04 
RMS  error  is:  0.90% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Partially 

Auto-covariance  -  Vertical  for  'Mwrun04' 

Mean  Standard  deviation  between  images  is:  11.38 
Power  Law  Model  Parameters: 
a_est  is:  115.44 
b_est  is:  0.00 
w_est  is:  0.03 
RMS  error  is:  1.10% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Partially 
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3.5.5  Image  Synthesis: 


Sample  Ensemble  Image  for  'Lwrun04' 
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Figure  25  -  Original  Images  from  Run  4 

Synthesized  Image  (Transformed)  for'Lwrun04' 
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Synthesized  Image  (Not  Transformed)  for  'Lwrun04' 
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Figure  26  -  Synthetic  Images  for  Run  4 


Sample  Ensemble  Image  for  'Mwrun04' 
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Synthesized  Image  (Transformed)  for  'Mwrun04' 
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Synthesized  Image  (Not  Transformed)  for  'Mwrun04' 
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3.6  Run  5 

3.6.1  General  Information 

Flight  Date:  30-6-97  NIGHT 
Altitude:  1.7K  ft 
Depression  Angle:  60  deg. 
Vegetation  Classes  in  Run:  11,4 
Predominate  Vegetation  Class:  54 

Ensemble  Size  is:  244 
Number  of  Valid  Images  is:  53 
Vertical  Resolution  is:  0.22m 
Horizontal  Resolution  is:  0.19m 


Histogram  Information  for  'Lwrun05' 

Normal  Model  Parameters: 

Mean  is:  1690.06 
Variance  is:  984.18 
RMS  error  is:  25.99% 

Chi2  P-value  is:  0.00% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 
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Histogram  Information  for  'Mwrun05' 

Normal  Model  Parameters: 

Mean  is:  683.10 
Variance  is:  281.89 
RMS  error  is:  42.06% 

Chi2  P-value  is:  0.00% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 


3.6.3  Auto-covariance  -  Horizontal 


Auto-covariance  -  Horizontal  Slice  for  'Lwnjn05‘  Auto-covariance  -  Horizontal  Slice  for  ’Mwrun05' 


Figure  28  -  Power  Law  Fit  -  Horizontal  Slice  of  Run  5 

Auto-covariance  -  Horizontal  for  'Lwrun05' 


Mean  Standard  deviation  between  images  is:  247.20 
Power  Law  Model  Parameters: 
a_est  is:  1627.02 
b_est  is:  1.04 
w_est  is:  0.03 
RMS  error  is:  1.13% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Partially 

Auto-covariance  -  Horizontal  for  'Mwrun05' 


Mean  Standard  deviation  between  images  is:  51.55 
Power  Law  Model  Parameters: 
a_est  is:  365.08 
b_est  is:  0.17 
w_est  is:  0.02 
RMS  error  is:  0.44% 

Model  Accepted  by  Threshold:  Yes 
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Model  Accepted  Visually:  Yes 

3.6.4  Auto-covariance  -  Vertical 


Auto-covariance  -  Vertical  Slice  for  'LwrunOS  Auto-covariance  -  Vertical  Slice  for  'Mwmn05' 


Figure  29  -  Power  Laiv  Fit  -  Vertical  Slice  of  Run  5 


Auto-covariance  -  Vertical  for  'Lwrun05' 

Mean  Standard  deviation  between  images  is:  252.79 
Power  Law  Model  Parameters: 
a_est  is:  1594.86 
b_est  is:  0.44 
w_est  is:  0.02 
RMS  error  is:  0.44% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Yes 

Auto-covariance  -  Vertical  for  'Mwrun05' 

Mean  Standard  deviation  between  images  is:  53.01 
Power  Law  Model  Parameters: 
a_est  is:  367.37 
b_est  is:  0.28 
w_est  is:  0.03 
RMS  error  is:  0.52% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Yes 
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3.6.5  Image  Synthesis: 


Sample  Ensemble  Image  for  'Lwrun05' 


Figure  30  -  Original  Images  from  Run  5 


Autocovariance  -  Vertical  Slice  for  ‘MwrunOS1 
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Synthesized  Image  (Transformed)  for  'Lwrun05' 
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Figure  31  -  Synthetic  Images  for  Run  5 


Synthesized  Image  (Transformed)  for'Mwrun05' 
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3.7  Run  6 

3.7.1  General  Information 

Flight  Date:  30-6-97  NIGHT 
Altitude:  1.7Kft 
Depression  Angle:  60  deg. 
Vegetation  Classes  in  Run:  48 
Predominate  Vegetation  Class:  48 

Ensemble  Size  is:  151 
Number  of  Valid  Images  is:  39 
Vertical  Resolution  is:  2.48m 
Horizontal  Resolution  is:  0.64m 


3.7.2  Histogram 


Histogram  for  'Lwrun06' 


Figure  32  -  Gaussian  Fit  -  Flistogram  of  Run  6 

Histogram  Information  for  'Lwrun06' 

Normal  Model  Parameters: 

Mean  is:  1683.04 
Variance  is:  234.60 
RMS  error  is:  9.94% 

Chi2  P-value  is:  0.00% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 


Histogram  for  'Mwrun06' 
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Histogram  Information  for  'Mwnm06' 

Normal  Model  Parameters: 

Mean  is:  688.01 
Variance  is:  70.00 
RMS  error  is:  18.25% 

Chi2  P-value  is:  0.00% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 


3.7.3  Auto-covariance  -  Horizontal 

Auto-covariance  -  Horizontal  Slice  for  'LwrunOS'  Auto-covariance  -  Horizontal  Slice  for  'MwrunOff 


Figure  33  -  Power  Law  Fit  -  Horizontal  Slice  of  Run  6 

Auto-covariance  -  Horizontal  for  'Lwrun06' 

Mean  Standard  deviation  between  images  is:  185.94 
Power  Law  Model  Parameters: 
a_est  is:  479.50 
b_est  is:  0.74 
w_est  is:  0.04 
RMS  error  is:  1.13% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  No 

Auto-covariance  -  Horizontal  for  'Mwrun06' 

Mean  Standard  deviation  between  images  is:  27.38 
Power  Law  Model  Parameters: 
a_est  is:  90.96 
b_est  is:  0.00 
w_est  is:  0.04 
RMS  error  is:  1.95% 

Model  Accepted  by  Threshold:  Yes 


34 


DSTO-TR-1456 


Model  Accepted  Visually:  No 

3.7.4  Auto-covariance  -  Vertical 

Auto-covariance  -  Vertical  Slice  for  ‘LwrunOff  Auto-covariance  -  Vertical  Slice  for  'Mwmn06' 


Figure  34  -  Power  Law  Fit  -  Vertical  Slice  of  Run  6 


Auto-covariance  -  Vertical  for  'Lwrun06' 

Mean  Standard  deviation  between  images  is:  184.95 
Power  Law  Model  Parameters: 
a_est  is:  452.59 
b_est  is:  0.00 
w_est  is:  0.03 
RMS  error  is:  1.35 
Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Partially 


Auto-covariance  -  Vertical  for  'Mwrun06' 

Mean  Standard  deviation  between  images  is:  27.53 
Power  Law  Model  Parameters: 
a_est  is:  89.54 
b_est  is:  0.00 
w_est  is:  0.04 
RMS  error  is:  1.37 
Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Partially 
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3.7.5  Image  Synthesis: 


Synthesized  Image  (Transformed)  for  'Lwrun06'  Synthesized  Image  (Transformed)  for  'Mwrun06' 
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Figure  36  -  Synthetic  Images  for  Run  6 
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3.8  Run  9 

3.8.1  General  Information 

Flight  Date:  1-7-97DAY 
Altitude:  lOKft 
Depression  Angle:  60  deg. 

Vegetation  Classes  in  Run:  105, 21, 11, 32,  4,  9 
Predominate  Vegetation  Class:  32 

Ensemble  Size  is:  150 
Number  of  Valid  Images  is:  47 
Vertical  Resolution  is:  1.30m 
Horizontal  Resolution  is:  1.13m 


3.8.2  Histogram 


Histogram  for  'Lwrun09' 


Figure  37  -  Gaussian  Fit  -  Histogram  of  Run 


Histogram  Information  for  'Lwrun09' 

Normal  Model  Parameters: 

Mean  is:  1963.37 
Variance  is:  14521.48 
RMS  error  is:  2.16% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  No 


Histogram  for  'Mwrun09' 
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Histogram  Information  for  'Mwrun09' 

Normal  Model  Parameters: 

Mean  is:  808.51 
Variance  is:  6007.26 
RMS  error  is:  3.70% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  No 

3.8.3  Auto-covariance  -  Horizontal 

Auto-covariance  -  Horizontal  Slice  for  'Lwrurr09'  Auto-covariance  -  Horizontal  Slice  for  'Mwrun09‘ 


Displacement  (m)  Displacement  (m) 


Figure  38  -  Poxver  Law  Fit  -  Horizontal  Slice  of  Run  9 

Auto-covariance  -  Horizontal  for  'Lwrun09' 

Mean  Standard  deviation  between  images  is:  907.60 
Power  Law  Model  Parameters: 
a_est  is:  14512.17 
b_est  is:  0.00 
w_est  is:  0.14 
RMS  error  is:  10.82% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 

Auto-covariance  -  Horizontal  for  'Mwrun09' 

Mean  Standard  deviation  between  images  is:  551.53 
Power  Law  Model  Parameters: 
a_est  is:  6714.88 
b  est  is:  0.00 
w_est  is:  0.11 
RMS  error  is:  5.48% 
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Auto-covariance  -  Vertical  Slice  for  'Lwnjn09'  Auto-covariance  -  Vertical  Slice  for  'Mwrun09' 


Figure  39  -  Poiver  Law  Fit  -  Vertical  Slice  of  Run  9 


Auto-covariance  -  Vertical  for  'Lwrun09' 

Mean  Standard  deviation  between  images  is:  883.06 
Power  Law  Model  Parameters: 
a_est  is:  9734.98 
b_est  is:  0.00 
w_est  is:  0.10 
RMS  error  is:  7.90% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 


Auto-covariance  -  Vertical  for  'Mwrun09' 

Mean  Standard  deviation  between  images  is:  518.84 
Power  Law  Model  Parameters: 
a_est  is:  6142.74 
b_est  is:  0.00 
w_est  is:  0.10 
RMS  error  is:  4.62% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 
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3.8.5  Image  Synthesis: 

Sample  Ensemble  Image  for  'LwrunOO'  Sample  Ensemble  Image  for  'Mwrun09' 
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Figure  40  -  Original  Images  from  Run  9 


Synthesized  Image  (Transformed)  for  'Lwrun09' 
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Figure  41  -  Synthetic  Images  for  Run  9 


Synthesized  Image  (T ransformed)  for  'Mwrun09‘ 
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3.9  Run  10 

3.9.1  General  Information 

Flight  Date:  -7-97  DAY 
Altitude:  OK  ft 
Depression  Angle:  60  deg. 

Vegetation  Classes  in  Run:  48, 11, 21, 105 
Predominate  Vegetation  Class:  11 

Ensemble  Size  is:  146 
Number  of  Valid  Images  is:  35 
Vertical  Resolution  is:  1.30m 
Horizontal  Resolution  is:  1.13m 


3.9.2  Histogram 


Histogram  for  'Lwrunl  0’ 


Histogram  Information  for  'LwrunlO' 

Normal  Model  Parameters: 

Mean  is:  1961.99 
Variance  is:  10000.00 
RMS  error  is:  2.20% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  No 


Histogram  for  ’MwrunlO' 
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Histogram  Information  for  'MwrunlO’ 

Normal  Model  Parameters: 

Mean  is:  812.46 
Variance  is:  3907.81 
RMS  error  is:  5.24% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  No 


3.9.3  Auto-covariance  -  Horizontal 


Auto-covariance  -  Horizontal  Slice  for  'LwrunlO'  Auto-covariance  -  Horizontal  Slice  for  'MwrunlO' 


Figure  43  -  Power  Law  Fit  -  Horizontal  Slice  of  Run  10 

Auto-covariance  -  Horizontal  for  'LwrunlO' 

Mean  Standard  deviation  between  images  is:  496.04 
Power  Law  Model  Parameters: 
a_est  is:  6591.87 
b_est  is:  0.00 
w_est  is:  0.10 
RMS  error  is:  6.27% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 

Auto-covariance  -  Horizontal  for  'MwrunlO' 

Mean  Standard  deviation  between  images  is:  305.81 
Power  Law  Model  Parameters: 
a_est  is:  3958.77 
b_est  is:  0.74 
w_est  is:  0.10 
RMS  error  is:  7.15% 
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Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 
3.9.4  Auto-covariance  -  Vertical 


Auto-covariance  -  Vertical  Slice  for  'Iwrunl O'  Auto-covariance  -  Vertical  Slice  for  'Mwrunl  O' 


Figure  44  -  Poxver  Laxv  Fit  -  Vertical  Slice  of  Run  10 


Auto-covariance  -  Vertical  for  'LwrunlO' 

Mean  Standard  deviation  between  images  is:  555.61 
Power  Law  Model  Parameters: 
a_est  is:  8157.81 
b_est  is:  0.00 
w_est  is:  0.12 
RMS  error  is:  6.12% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 


Auto-covariance  -  Vertical  for  'MwrunlO' 

Mean  Standard  deviation  between  images  is:  288.02 
Power  Law  Model  Parameters: 
a_est  is:  4670.79 
b_est  is:  1.25 
w_est  is:  0.13 
RMS  error  is:  5.01% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Yes 
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3.9.5  Image  Synthesis: 


Sample  Ensemble  Image  for  'LwrunlO'  Sample  Ensemble  Image  for  'MwrunlO' 
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Figure  45  -  Original  Images  from  Run  10 


Synthesized  Image  (Transformed)  for  'LwrunlO' 
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Figure  46  -  Synthetic  Inmgesfor  Run  10 
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3.10  Run  11 

3.10.1  General  Information 

Flight  Date:  3-7-97  DAY 
Altitude:  10K  ft 
Depression  Angle:  60  deg. 
Vegetation  Classes  in  Run:  21, 15, 51 
Predominate  Vegetation  Class:  21 

Ensemble  Size  is:  150 
Number  of  Valid  Images  is:  17 
Vertical  Resolution  is:  1.30m 
Horizontal  Resolution  is:  1.13m 


3.10.2  Histogram 


Histogram  for  'Lwrunl  1 '  x  1 0‘3  Histogram  for  'Mwrunl  1 ' 


Figure  47  -  Gaussian  Fit  -  Histogram  of  Run  11 


Histogram  Information  for  'Lwrunll' 

Normal  Model  Parameters: 

Mean  is:  2118.03 
Variance  is:  27839.50 
RMS  error  is:  2.16% 

Chi2P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  No 
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Histogram  Information  for  'Mwrunll' 

Normal  Model  Parameters: 

Mean  is:  917.71 
Variance  is:  12503.97 
RMS  error  is:  3.14% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  No 

3.10.3  Auto-covariance  -  Horizontal 

v  m4  Auto-covariance  -  Horizontal  Slice  for  'LwrunIT  Auto-covariance  -  Horizontal  Slice  for  'Mwrunl  V 


Figure  48  -  Power  Law  Fit  -  Horizontal  Slice  of  Run  11 

Auto-covariance  -  Horizontal  for  'Lwrunll' 

Mean  Standard  deviation  between  images  is:  2527.34 
Power  Law  Model  Parameters: 
a_est  is:  18431.15 
b_est  is:  0.00 
w_est  is:  0.11 
RMS  error  is:  7.50% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 

Auto-covariance  -  Horizontal  for  'Mwrunll' 

Mean  Standard  deviation  between  images  is:  1080.30 
Power  Law  Model  Parameters: 
a_est  is:  9711.80 
b_est  is:  0.00 
w_est  is:  0.12 
RMS  error  is:  9.10% 
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Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 
3.10.4  Auto-covariance  -  Vertical 
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x  Auto-covariance  -  Vertical  Slice  for  twrunl  1 '  Auto-covariance  -  Vertical  Slice  for  'Mwrunl  1 ' 


Figure  49  -  Poiver  Law  Fit  -  Vertical  Slice  of  Run  11 

Auto-covariance  -  Vertical  for  'Lwrunll' 

Mean  Standard  deviation  between  images  is:  2616.55 
Power  Law  Model  Parameters: 
a_est  is:  18068.13 
b_est  is:  0.00 
w_est  is:  0.12 
RMS  error  is:  5.70% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 


Auto-covariance  -  Vertical  for  'Mwrunll' 

Mean  Standard  deviation  between  images  is:  1145.11 
Power  Law  Model  Parameters: 
a_est  is:  9941.85 
b_est  is:  0.00 
w_est  is:  0.13 
RMS  error  is:  7.07% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 
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3.10.5  Image  Synthesis: 


Sample  Ensemble  Image  for  'Lwrunl  1’ 


Sample  Ensemble  Image  for  'Mwrunli' 


Figure  50  -  Original  Images  from  Run  11 


Synthesized  Image  (Transformed)  for 'Lwrunl  1' 
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Figure  51  -  Synthetic  Images  for  Run  11 


Synthesized  Image  (Transformed)  for  'Mwrunli' 


20 

40 

60 

80 

100 

120 


50  100  150  200  250 


48 


DSTO-TR-1456 


3.11  Run  12 

3.11.1  General  Information 

Flight  Date:  3-7-97DAY 
Altitude:  lOKft 
Depression  Angle:  60  deg. 

Vegetation  Classes  in  Run:  13,  51, 15, 4, 21 
Predominate  Vegetation  Class:  15 

Ensemble  Size  is:  269 
Number  of  Valid  Images  is:  31 
Vertical  Resolution  is:  1.30m 
Horizontal  Resolution  is:  1.13m 

3.11.2  Histogram 

x-]0'3  Histogram  for 'Lwrunl 2'  xIO"3  Histogram  for 'Mwrun12' 


Figure  52  -  Gaussian  Fit  -  Histogram  of  Run  12 

Histogram  Information  for  'Lwrunl2' 

Normal  Model  Parameters: 

Mean  is:  2081.22 
Variance  is:  22712.26 
RMS  error  is:  1.29% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  No 
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Histogram  Information  for  'Mwrunl2' 

Normal  Model  Parameters: 

Mean  is:  875.47 
Variance  is:  9580.54 
RMS  error  is:  2.33% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Partially 

3.11.3  Auto-covariance  -  Horizontal 


/Auto-covariance  -  Horizontal  Slice  for  'Lwrunf  2'  Auto-covariance  -  Horizontal  Slice  for  'Mwrun12' 


Figure  53  -  Power  Laxo  Fit  -  Horizontal  Slice  of  Run  12 

Auto-covariance  -  Horizontal  for  'Lwrunl2' 

Mean  Standard  deviation  between  images  is:  1369.61 
Power  Law  Model  Parameters: 
a_est  is:  11251.79 
b_est  is:  0.00 
w_est  is:  0.08 
RMS  error  is:  7.96% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 

Auto-covariance  -  Horizontal  for  'Mwrunl2' 

Mean  Standard  deviation  between  images  is:  610.68 
Power  Law  Model  Parameters: 
a_est  is:  6331.00 
b_est  is:  0.00 
w_est  is:  0.09 
RMS  error  is:  9.17% 
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Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 
3.11.4  Auto-covariance  -  Vertical 
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Auto-covariance  -  Vertical  Slice  for  'Mwrunl  2' 


Figure  54  -  Power  Law  Fit  -  Vertical  Slice  of  Run  12 


Auto-covariance  -  Vertical  for  'Lwrunl2' 

Mean  Standard  deviation  between  images  is:  1348.69 
Power  Law  Model  Parameters: 
a_est  is:  11000.20 
b_est  is:  0.00 
w_est  is:  0.09 
RMS  error  is:  8.77% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 

Auto-covariance  -  Vertical  for  'Mwrunl2' 

Mean  Standard  deviation  between  images  is:  598.14 
Power  Law  Model  Parameters: 
a_est  is:  6469.22 
b_est  is:  0.00 
w_est  is:  0.11 
RMS  error  is:  6.85% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 
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3.11.5  Image  Synthesis: 
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Figure  55  -  Original  Images  from  Run  12 


Synthesized  Image  (Transformed)  for'Lwrun12' 
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Figure  56  -  Synthetic  Images  for  Run  12 


Synthesized  Image  (Transformed)  for 'Mwrunl 2' 
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3.12  Run  13 

3.12.1  General  Information 

Flight  Date:  3-7-97  DAY 
Altitude:  5Kft 
Depression  Angle:  60  deg. 

Vegetation  Classes  in  Run:  9, 21, 11, 32 
Predominate  Vegetation  Class:  32 

Ensemble  Size  is:  250 
Number  of  Valid  Images  is:  62 
Vertical  Resolution  is:  0.65m 
Horizontal  Resolution  is:  0.56m 


3.12.2  Histogram 


Histogram  for  'Lwrunl  3' 


Greylevel 

Figure  57  -  Gaussian  Fit  -  Histogram  of  Run  13 


Histogram  for  ’Mwrun13' 


Histogram  Information  for  'Lwrunl3' 

Normal  Model  Parameters: 

Mean  is:  2066.47 
Variance  is:  21795.67 
RMS  error  is:  3.45% 

Chi2  P-value  is:  99.55% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  No 
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Histogram  Information  for  'Mwrunl3' 

Normal  Model  Parameters: 

Mean  is:  904.15 
Variance  is:  10624.08 
RMS  error  is:  4.31% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  No 

3.12.3  Auto-covariance  -  Horizontal 

x  iq4  Auto-covariance  -  Horizontal  Slice  for  'Lwrunl 3'  AutMovariance  -  Horizontal  Slice  for  'Mwnin13' 


Figure  58  -  Power  Laxv  Fit  -  Horizontal  Slice  of  Run  13 

Auto-covariance  -  Horizontal  for  'Lwrunl3' 

Mean  Standard  deviation  between  images  is:  666.45 
Power  Law  Model  Parameters: 
a_est  is:  8080.16 
b_est  is:  0.00 
w_est  is:  0.09 
RMS  error  is:  15.29% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 

Auto-covariance  -  Horizontal  for  'Mwrunl3' 

Mean  Standard  deviation  between  images  is:  392.46 
Power  Law  Model  Parameters: 
a_est  is:  8413.28 
b_est  is:  0.00 
w_est  is:  0.15 
RMS  error  is:  14.47% 
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Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 
3.12.4  Auto-covariance  -  Vertical 
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x  iq4  Auto-covariance  -  Vertical  Slice  for  'Lwnjn13'  Auto-covariance  -  Vertical  Slice  for  'Mwrun13' 


Figure  59  -  Power  Law  Fit  -  Vertical  Slice  of  Run  13 


Auto-covariance  -  Vertical  for  'Lwrunl3' 

Mean  Standard  deviation  between  images  is:  641.98 
Power  Law  Model  Parameters: 
a_est  is:  7837.69 
b_est  is:  0.00 
w_est  is:  0.10 
RMS  error  is:  20.77% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 

Auto-covariance  -  Vertical  for  'Mwrunl3' 

Mean  Standard  deviation  between  images  is:  390.10 
Power  Law  Model  Parameters: 
a_est  is:  7605.90 
b_est  is:  0.00 
w_est  is:  0.16 
RMS  error  is:  16.36% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 
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3.12.5  Image  Synthesis: 

Sample  Ensemble  lmagefor'Lwrun13'  Sample  Ensemble  Image  for  'Mwrun13' 


Figure  60  -  Original  Images  from  Run  13 


Synthesized  Image  (Transformed)  for  'Lwrun13’ 


Figure  61  -  Synthetic  Images  for  Run  13 


Synthesized  Image  (Transformed)  for  'Mwrun13' 
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3.13  Run  14 

3.13.1  General  Information 

Flight  Date:  3-7-97DAY 
Altitude:  5Kft 
Depression  Angle:  60  deg. 
Vegetation  Classes  in  Run:  48, 11, 51 
Predominate  Vegetation  Class  11 

Ensemble  Size  is:  250 
Number  of  Valid  Images  is:  75 
Vertical  Resolution  is:  0.65m 
Horizontal  Resolution  is:  0.56m 


3.13.2  Histogram 


1800  2000  2200  2400  2600  800  1000  1200  1400 

Greylevel  Grey  level 

Figure  62  -  Gaussian  Fit  -  Histogram  of  Run  14 

Histogram  Information  for  'Lwrunl4' 

Normal  Model  Parameters: 

Mean  is:  2147.89 
Variance  is:  24725.98 
RMS  error  is:  2.26% 

Chi2  P- value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  No 
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Histogram  Information  for  'Mwrunl4' 

Normal  Model  Parameters: 

Mean  is:  967.89 
Variance  is:  10878.77 
RMS  error  is:  2.86% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  No 

3.13.3  Auto-covariance  -  Horizontal 


Displacement  (m)  Displacement  (m) 

Figure  63  -  Power  Law  Fit  -  Horizontal  Slice  of  Run  14 

Auto-covariance  -  Horizontal  for  'Lwrunl4' 

Mean  Standard  deviation  between  images  is:  856.52 
Power  Law  Model  Parameters: 
a_est  is:  22984.96 
b_est  is:  12.00 
w_est  is:  0.20 
RMS  error  is:  12.57% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 

Auto-covariance  -  Horizontal  for  'Mwrunl4' 

Mean  Standard  deviation  between  images  is:  473.18 
Power  Law  Model  Parameters: 
a_est  is:  8896.77 
b_est  is:  0.19 
w_est  is:  0.17 
RMS  error  is:  6.44% 
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Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Yes 
3.13.4  Auto-covariance  -  Vertical 
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Auto-covariance  -  Vertical  Slice  for  'Mwrun14' 


Figure  64  -  Poxver  Law  Fit  -  Vertical  Slice  of  Run  14 


Auto-covariance  -  Vertical  for  'Lwrunl4' 

Mean  Standard  deviation  between  images  is:  784.67 
Power  Law  Model  Parameters: 
a_est  is:  24578.28 
b_est  is:  0.68 
w_est  is:  0.23 
RMS  error  is:  4.89% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Yes 

Auto-covariance  -  Vertical  for  'Mwrunl4' 

Mean  Standard  deviation  between  images  is:  441.07 
Power  Law  Model  Parameters: 
a_est  is:  16219.27 
b_est  is:  2.07 
w_est  is:  0.26 
RMS  error  is:  5.68% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Yes 
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3.13.5  Image  Synthesis: 


Sample  Ensemble  Image  for  'Lwrun14'  Sample  Ensemble  Image  for  'Mwrun14' 


50  100  150  200  250  50  100  150  200  250 


Figure  65  -  Original  Images  from  Run  14 


Synthesized  Image  (Transformed)  for  ‘Mwrun14‘ 


50  100  150  200  250 


Figure  66  -  Synthetic  Images  for  Run  14 


60 


DSTO-TR-1456 


3.14  Run  15 

3.14.1  General  Information 

Flight  Date:  3-7-97DAY 
Altitude:  5K  ft 
Depression  Angle:  60  deg. 
Vegetation  Classes  in  Run:  21, 15, 51 
Predominate  Vegetation  Class:  21 

Ensemble  Size  is:  250 
Number  of  Valid  Images  is:  49 
Vertical  Resolution  is:  0.65m 
Horizontal  Resolution  is:  0.56m 


3.14.2  Histogram 


Histogram  for  'Lwrun15' 


Figure  67  -  Gaussian  Fit  -  Histogram  of  Run  15 


Histogram  for  'Mwmn15' 


Greylevel 


Histogram  Information  for  'Lwrunl5' 

Normal  Model  Parameters: 

Mean  is:  2198.06 
Variance  is:  28755.24 
RMS  error  is:  3.95% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  No 
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Histogram  Information  for  'Mwrunl5' 

Normal  Model  Parameters: 

Mean  is:  1044.27 
Variance  is:  13828.69 
RMS  error  is:  6.07% 

Chi2  P-value  is:  0.02% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 

3.14.3  Auto-covariance  -  Horizontal 

x  10<  Auto-covariance  -  Horizontal  Slice  for  'Lwrun15'  x  i04Autooovariance  -  Horizontal  Slice  for  'Mwrunl  5' 


Figure  68  -  Power  Laxv  Fit  -  Horizontal  Slice  of  Run  15 

Auto-covariance  -  Horizontal  for  'Lwrunl5' 

Mean  Standard  deviation  between  images  is:  1523.55 
Power  Law  Model  Parameters: 
a_est  is:  28049.33 
b_est  is:  0.00 
w_est  is:  0.12 
RMS  error  is:  6.68% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 

Auto-covariance  -  Horizontal  for  'Mwrunl5' 

Mean  Standard  deviation  between  images  is:  880.74 
Power  Law  Model  Parameters: 
a_est  is:  16413.18 
b_est  is:  0.00 
w_est  is:  0.12 
RMS  error  is:  5.11% 
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Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 
3.14.4  Auto-covariance  -  Vertical 
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x  i  Q4  Auto-covariance  -  Vertical  Slice  for  ‘Lwrunl  5' 


Figure  69  -  Power  Law  Fit  -  Vertical  Slice  of  Run  15 


Auto-covariance  -  Vertical  for  'Lwrunl5' 

Mean  Standard  deviation  between  images  is:  1431.44 
Power  Law  Model  Parameters: 
a_est  is:  25309.82 
b_est  is:  0.00 
w_est  is:  0.11 
RMS  error  is:  5.64% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 

Auto-covariance  -  Vertical  for  'Mwrunl5' 

Mean  Standard  deviation  between  images  is:  841.12 
Power  Law  Model  Parameters: 
a_est  is:  15107.65 
b_est  is:  0.00 
w_est  is:  0.12 
RMS  error  is:  4.88% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 
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3.14.5  Image  Synthesis: 


Sample  Ensemble  Image  for 'Lwrun15'  Sample  Ensemble  Image  for  'MwrunlS' 


Figure  70  -  Original  Images  from  Run  15 


Synthesized  Image  (Transformed)  for  'Lwrunl 5'  Synthesized  Image  (Transformed)  for  'Mwrunl 5‘ 


Figure  71  -  Synthetic  Images  for  Run  15 
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3.15  Run  16 

3.15.1  General  Information 

Flight  Date:  3-7-97DAY 
Altitude:  5K  ft 
Depression  Angle:  60  deg. 

Vegetation  Classes  in  Run:  13,  51, 15, 4,  21 
Predominate  Vegetation  Class:  15 

Ensemble  Size  is:  548 
Number  of  Valid  Images  is:  60 
Vertical  Resolution  is:  0.65m 
Horizontal  Resolution  is:  0.56m 


3.15.2  Histogram 


Histogram  for  'Lwrunl  6' 


Figure  72  -  Gaussian  Fit  -  Histogram  of  Run  16 


3  Histogram  for  'Mwrunl  6' 


Greylevel 


Histogram  Information  for  'Lwrunl6' 

Normal  Model  Parameters: 

Mean  is:  2182.96 
Variance  is:  34618.21 
RMS  error  is:  1.63% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  No 
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Histogram  Information  for  'Mwrunl6' 

Normal  Model  Parameters: 

Mean  is:  1007.99 
Variance  is:  16445.23 
RMS  error  is:  2.20% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  No 


3.15.3  Auto-covariance  -  Horizontal 


Auto-covariance  -  Horizontal  Slice  for  'Mwrun16' 


Displacement  (m) 


Figure  73  -  Power  Law  Fit  -  Horizontal  Slice  of  Run  16 


Auto-covariance  -  Horizontal  for  'Lwrunl6' 

Mean  Standard  deviation  between  images  is:  1242.94 
Power  Law  Model  Parameters: 
a_est  is:  13110.58 
b_est  is:  0.00 
w_est  is:  0.11 
RMS  error  is:  12.09% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 

Auto-covariance  -  Horizontal  for  'Mwrunl6' 

Mean  Standard  deviation  between  images  is:  654.40 
Power  Law  Model  Parameters: 
a_est  is:  8054.12 
b_est  is:  0.00 
w_est  is:  0.12 
RMS  error  is:  10.14% 
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Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 
3.15.4  Auto-covariance  -  Vertical 
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Displacement  (m)  Displacement  (m) 


Figure  74  -  Poiver  Law  Fit  -  Vertical  Slice  of  Run  16 

Auto-covariance  -  Vertical  for  'Lwrunl6' 

Mean  Standard  deviation  between  images  is:  1197.24 
Power  Law  Model  Parameters: 
a_est  is:  16968.82 
b_est  is:  0.00 
w_est  is:  0.15 
RMS  error  is:  10.49% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 

Auto-covariance  -  Vertical  for  'Mwrunl6' 

Mean  Standard  deviation  between  images  is:  630.76 
Power  Law  Model  Parameters: 
a_est  is:  10116.72 
b_est  is:  0.00 
w_est  is:  0.16 
RMS  error  is:  8.24% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 
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3.15.5  Image  Synthesis: 


Sample  Ensemble  Image  for ’Lwrun16‘ 
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Figure  75  -  Original  Images  from  Run  16 


Sample  Ensemble  Image  for  'Mwrun16’ 
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Synthesized  Image  (Transformed)  for  'Lwrunl  6' 
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Figure  76  -  Synthetic  Images  for  Run  16 


Synthesized  Image  (Transformed)  for  'Mwrun16' 
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3.16  Run  17 

3.16.1  General  Information 

Flight  Date:  3-7-97  DAY 
Altitude:  10K  ft 
Depression  Angle:  60  deg. 
Vegetation  Classes  in  Run:  4, 106 
Predominate  Vegetation  Class:  106 

Ensemble  Size  is:  200 
Number  of  Valid  Images  is:  50 
Vertical  Resolution  is:  1.30m 
Horizontal  Resolution  is:  1.13m 


3.16.2  Histogram 


Histogram  for  'Lwrun17' 


Figure  77  -  Gaussian  Fit  -  Histogram  of  Run 


Histogram  for  'Mwrun17' 


Grey  level 


Histogram  Information  for  'Lwrunl7' 

Normal  Model  Parameters: 

Mean  is:  1898.02 
Variance  is:  4800.00 
RMS  error  is:  2.24% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Yes 
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Histogram  Information  for  'Mwrunl7' 

Normal  Model  Parameters: 

Mean  is:  735.49 
Variance  is:  1473.30 
RMS  error  is:  4.73% 

Chi2  P-value  is:  0.00% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 

3.16.3  Auto-covariance  -  Horizontal 


Auto-covariance  -  Horizontal  Slice  for  'Lwrunl  T  Auto-covariance  -  Horizontal  Slice  for  'Mwrunl  T 


Figure  78  -  Power  Law  Fit  -  Horizontal  Slice  of  Run  17 

Auto-covariance  -  Horizontal  for  'Lwrunl7' 

Mean  Standard  deviation  between  images  is:  546.37 
Power  Law  Model  Parameters: 
a__est  is:  5724.68 
b_est  is:  0.00 
w_estis:  0.11 
RMS  error  is:  6.50% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 

Auto-covariance  -  Horizontal  for  'Mwrunl7' 

Mean  Standard  deviation  between  images  is:  251.93 
Power  Law  Model  Parameters: 
a_est  is:  2248.62 
b_est  is:  0.00 
w_est  is:  0.11 
RMS  error  is:  7.60% 
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Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 
3.16.4  Auto-covariance  -  Vertical 


DSTO-TR-1456 


Auto-covariance  -  Vertical  Slice  for  'Lwrun17' 


Auto-covariance  -  Vertical  Slice  for  ’Mwrunl  T 


Figure  79  -  Power  Law  Fit  -  Vertical  Slice  of  Run  17 


Auto-covariance  -  Vertical  for  'Lwrunl7' 

Mean  Standard  deviation  between  images  is:  560.96 
Power  Law  Model  Parameters: 
a_est  is:  4437.64 
b_est  is:  0.00 
w_est  is:  0.08 
RMS  error  is:  4.95% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 

Auto-covariance  -  Vertical  for  'Mwrunl7' 

Mean  Standard  deviation  between  images  is:  255.09 
Power  Law  Model  Parameters: 
a_est  is:  2223.63 
b_est  is:  0.00 
w_est  is:  0.11 
RMS  error  is:  3.99% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 
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3.16.5  Image  Synthesis: 


Sample  Ensemble  Image  for  'Lwrun17'  Sample  Ensemble  Image  for  'Mwrun17' 
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Figure  80  -  Original  Images  from  Run  1 7 


Figure  81  -  Synthetic  Images  for  Run  17 
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3.17  Run  18 


3.17.1  General  Information 

Flight  Date:  3-7-97  DAY 
Altitude:  10K  ft 
Depression  Angle:  60  deg. 
Vegetation  Classes  in  Run:  53, 4, 106 
Predominate  Vegetation  Class:  53 

Ensemble  Size  is:  251 
Number  of  Valid  Images  is:  26 
Vertical  Resolution  is:  1.30m 
Horizontal  Resolution  is:  1.13m 


3.17.2  Histogram 


Histogram  for  'Lwrunl  8' 


Greylevel 


Figure  82  -  Gaussian  Fit  -  Histogram  of  Run  18 

Histogram  Information  for  'Lwrunl8' 

Normal  Model  Parameters: 

Mean  is:  1931.10 
Variance  is:  4800.00 
RMS  error  is:  4.56% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  No 


Histogram  for  'MwrunIS' 
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Histogram  Information  for  ’Mwrunl8' 

Normal  Model  Parameters: 

Mean  is:  744.81 
Variance  is:  1663.81 
RMS  error  is:  6.01  % 

CM2  P-value  is:  0.13% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 

3.17.3  Auto-covariance  -  Horizontal 

Auto-covariance  -  Horizontal  Slice  for  'Lwrunl  8'  Auto-covariance  -  Horizontal  Slice  for  'Mwrunl  8' 


Figure  83  -  Poioer  Law  Fit  -  Horizontal  Slice  of  Run  18 

Auto-covariance  -  Horizontal  for  'Lwrunl8' 

Mean  Standard  deviation  between  images  is:  1036.54 
Power  Law  Model  Parameters: 
a_est  is:  7730.06 
b_est  is:  0.00 
w_est  is:  0.18 
RMS  error  is:  13.85% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 

Auto-covariance  -  Horizontal  for  'Mwrunl8' 

Mean  Standard  deviation  between  images  is:  424.45 
Power  Law  Model  Parameters: 
a_est  is:  3631.64 
b_est  is:  2.40 
w_est  is:  0.18 
RMS  error  is:  16.44% 

Model  Accepted  by  Threshold:  No 
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Model  Accepted  Visually:  Partially 

3.17.4  Auto-covariance  -  Vertical 


Auto-covariance  -  Vertical  Slice  for  'Lv/run18'  Auto-covariance  -  Vertical  Slice  for  'Mwnjn18' 


Figure  84  -  Power  Law  Fit  -  Vertical  Slice  of  Run  18 


Auto-covariance  -  Vertical  for  'Lwrunl8' 

Mean  Standard  deviation  between  images  is:  1007.21 
Power  Law  Model  Parameters: 
a_est  is:  6164.67 
b_est  is:  0.00 
w_est  is:  0.13 
RMS  error  is:  5.47% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Yes 


Auto-covariance  -  Vertical  for  'Mwrunl8' 

Mean  Standard  deviation  between  images  is:  406.39 
Power  Law  Model  Parameters: 
a_est  is:  3113.55 
b_est  is:  1.67 
w_est  is:  0.14 
RMS  error  is:  5.00% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Yes 
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3.17.5  Image  Synthesis: 


Sample  Ensemble  Image  for  ‘Lwrunl  8'  Sample  Ensemble  li 
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Figure  85  -  Original  Images  from  Run  18 


Synthesized  Image  (Transformed)  for  'Lwrun18' 
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Figure  86  -  Synthetic  Images  for  Run  18 


Synthesized  Image  (Transformed)  for  'Mwrun18' 
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3.18  Run  19 

3.18.1  General  Information 

Flight  Date:  3-7-97DAY 
Altitude:  5Kft 
Depression  Angle:  60  deg. 
Vegetation  Classes  in  Run:  4, 106 
Predominate  Vegetation  Class:  106 

Ensemble  Size  is:  414 
Number  of  Valid  Images  is:  103 
Vertical  Resolution  is:  0.65m 
Florizontal  Resolution  is:  0.56m 


3.18.2  Histogram 


Histogram  for  'Lwrunl  9'  Histogram  for  'Mwrun19' 


Histogram  Information  for  'Lwrunl  9' 

Normal  Model  Parameters: 

Mean  is:  1942.97 
Variance  is:  3500.00 
RMS  error  is:  4.03% 

Chi2  P-value  is:  89.65% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 


77 


DSTO-TR-1456 


Histogram  Information  for  'Mwmnl9' 

Normal  Model  Parameters: 

Mean  is:  823.38 
Variance  is:  1600.00 
RMS  error  is:  5.09% 

Chi2  P-value  is:  0.00% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 

3.18.3  Auto-covariance  -  Horizontal 

Auto-covariance  -  Horizontal  Slice  for  'Lwrun20'  Auto-covariance  -  Horizontal  Slice  for  'Mwrun20' 


Figure  88  -  Power  Law  Fit  -  Horizontal  Slice  of  Run  19 

Auto-covariance  -  Horizontal  for  'Lwrunl9' 

Mean  Standard  deviation  between  images  is:  186.04 
Power  Law  Model  Parameters: 
a_est  is:  4792.95 
b_est  is:  0.00 
w_est  is:  0.18 
RMS  error  is:  4.10% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Yes 

Auto-covariance  -  Horizontal  for  'Mwrunl9' 

Mean  Standard  deviation  between  images  is:  101.75 
Power  Law  Model  Parameters: 
a_est  is:  2483.05 
b_est  is:  0.00 
w_est  is:  0.17 
RMS  error  is:  3.54% 
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Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 
3.18.4  Auto-covariance  -  Vertical 
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Auto-covariance  -  Vertical  Slice  for  'Lwrun19'  Auto-covariance  -  Vertical  Slice  for  'Mwrun19' 


Figure  89  -  Power  Law  Fit  -  Vertical  Slice  of  Run  19 


Auto-covariance  -  Vertical  for  'Lwrunl9' 

Mean  Standard  deviation  between  images  is:  181.17 
Power  Law  Model  Parameters: 
a_est  is:  3352.59 
b_est  is:  0.00 
w_est  is:  0.14 
RMS  error  is:  7.81% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 

Auto-covariance  -  Vertical  for  'Mwrunl9' 

Mean  Standard  deviation  between  images  is:  101.26 
Power  Law  Model  Parameters: 
a_est  is:  1970.02 
b_est  is:  0.00 
w_est  is:  0.14 
RMS  error  is:  4.61% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 
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3.18.5  Image  Synthesis: 


Sample  Ensemble  Image  for'Lwrun19' 
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Figure  90  -  Original  linages  from  Run  19 


Synthesized  Image  (Transformed)  for'Lwrun19' 
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Figure  91  -  Synthetic  Images  for  Run  19 
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Synthesized  Image  (Transformed)  for  'Mwrun19' 
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3.19  Run  20 

3.19.1  General  Information 

Flight  Date:  3-7-97DAY 
Altitude:  5K  ft 
Depression  Angle:  60  deg. 
Vegetation  Classes  in  Run:  53,  4, 106 
Predominate  Vegetation  Class:  53 

Ensemble  Size  is:  495 
Number  of  Valid  Images  is:  65 
Vertical  Resolution  is:  0.65m 
Horizontal  Resolution  is:  0.56m 


3.19.2  Histogram 


Histogram  for  'Lwrun20'  Histogram  for  'Mwrun20' 


Histogram  Information  for  'Lwrun20' 

Normal  Model  Parameters: 

Mean  is:  1947.65 
Variance  is:  3100.00 
RMS  error  is:  5.28% 

Chi2  P-value  is:  0.00% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 
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Histogram  Information  for  'Mwrun20' 

Normal  Model  Parameters: 

Mean  is:  817.25 
Variance  is:  1350.00 
RMS  error  is:  8.38% 

Chi2  P-value  is:  0.00% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 

3.19.3  Auto-covariance  -  Horizontal 

Auto-covariance  -  Horizontal  Slice  For  'Lwwn20’  Autcvcovariance  -  Horizontal  Slice  for  'Mwrun20' 


Figure  93  -  Poiver  Law  Fit  -  Horizontal  Slice  of  Run  20 

Auto-covariance  -  Horizontal  for  'Lwrun20' 

Mean  Standard  deviation  between  images  is:  211.92 
Power  Law  Model  Parameters: 
a_est  is:  6632.03 
b_est  is:  0.00 
w_est  is:  0.25 
RMS  error  is:  9.38% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Yes 

Auto-covariance  -  Horizontal  for  'Mwrun20' 

Mean  Standard  deviation  between  images  is:  131.21 
Power  Law  Model  Parameters: 
a_est  is:  2656.06 
b_est  is:  0.00 
w_est  is:  0.19 
RMS  error  is:  12.69% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 
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3.19.4  Auto-covariance  -  Vertical 

Auto-covariance  -  Vertical  Slice  for  'Lwnjn20'  Auto-covariance  -  Vertical  Slice  for  'Mwrun20' 


Figure  94  -  Power  Law  Fit  -  Vertical  Slice  of  Run  20 


Auto-covariance  -  Vertical  for  ’Lwrun20' 

Mean  Standard  deviation  between  images  is:  214.62 
Power  Law  Model  Parameters: 
a_est  is:  5603.16 
b_est  is:  0.00 
w_est  is:  0.22 
RMS  error  is:  11.43% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Yes 

Auto-covariance  -  Vertical  for  'Mwrun20' 

Mean  Standard  deviation  between  images  is:  134.95 
Power  Law  Model  Parameters: 
a_est  is:  2612.09 
b_est  is:  0.00 
w_est  is:  0.19 
RMS  error  is:  12.77% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 
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3.19.5  Image  Synthesis: 


Sample  Ensemble  Image  for  'Lwrun20' 
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Figure  95  -  Original  Images  from  Run  20 


Sample  Ensemble  Image  for  ’Mwrun20' 
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Synthesized  Image  (Transformed)  for  'Lwrijn20’  Synthesized  Image  (Transformed)  for  'Mwrun20' 
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Figure  96  -  Syntltetic  Images  for  Run  20 
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3.20  Run  21 

3.20.1  General  Information 

Flight  Date:  3-7-97DAY 
Altitude:  lOKft 
Depression  Angle:  60  deg. 
Vegetation  Classes  in  Run:  18, 106 
Predominate  Vegetation  Class:  106 

Ensemble  Size  is:  150 
Number  of  Valid  Images  is:  73 
Vertical  Resolution  is:  1.30m 
Horizontal  Resolution  is:  1.13m 


3.20.2  Histogram 

Histogram  for  'Lwrun21'  Histogram  for  'Mwrun21‘ 


Figure  97  -  Gaussian  Fit  -  Flistogram  of  Run  21 


Histogram  Information  for  'Lwrun21' 

Normal  Model  Parameters: 

Mean  is:  1738.81 
Variance  is:  680.00 
RMS  error  is:  6.88% 

Chi2  P-value  is:  0.00% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 
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Histogram  Information  for  'Mwrun21' 

Normal  Model  Parameters: 

Mean  is:  632.47 
Variance  is:  201.06 
RMS  error  is:  15.70% 

Chi2  P-value  is:  0.00% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 

3.20.3  Auto-covariance  -  Horizontal 


Auto-covariance  -  Horizontal  Slice  for  'Lwrun21 '  Auto-covariance  -  Horizontal  Slice  for  'Mwrun21 ' 


Figure  98  -  Power  Laxo  Fit  -  Horizontal  Slice  of  Run  21 

Auto-covariance  -  Horizontal  for  'Lwrun21' 

Mean  Standard  deviation  between  images  is:  94.50 
Power  Law  Model  Parameters: 
a_est  is:  951.44 
b_est  is:  0.00 
w_est  is:  0.08 
RMS  error  is:  3.13% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Yes 

Auto-covariance  -  Horizontal  for  'Mwrun21' 

Mean  Standard  deviation  between  images  is:  21.47 
Power  Law  Model  Parameters: 
a_est  is:  274.83 
b_est  is:  0.00 
w_est  is:  0.08 
RMS  error  is:  3.10% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Yes 
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3.20.4  Auto-covariance  -  Vertical 


Auto-covariance  -  Horizontal  Slice  for  ’Lwrun21' 


Auto-covariance  -  Horizontal  Slice  for  'Mwrun21' 


Figure  99  -  Power  Laiv  Fit  -  Vertical  Slice  of  Run  21 

Auto-covariance  -  Vertical  for  'Lwrun21' 

Mean  Standard  deviation  between  images  is:  89.90 
Power  Law  Model  Parameters: 
a_est  is:  971.88 
b_est  is:  0.00 
w_est  is:  0.09 
RMS  error  is:  3.27% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Yes 


Auto-covariance  -  Vertical  for  ’Mwrun21' 

Mean  Standard  deviation  between  images  is:  20.34 
Power  Law  Model  Parameters: 
a_est  is:  275.17 
b_est  is:  0.00 
w_est  is:  0.09 
RMS  error  is:  2.33% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Yes 
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3.20.5  Image  Synthesis: 
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3.21  Runs  22-23 

3.21.1  General  Information 

Flight  Date:  3-7-97  DAY 
Altitude:  10K  ft 
Depression  Angle:  60  deg. 

Vegetation  Classes  in  Run:  106, 3, 47,  3, 105 
Predominate  Vegetation  Class:  3 

Ensemble  Size  is:  233 
Number  of  Valid  Images  is:  233 
Vertical  Resolution  is:  1.30m 
Horizontal  Resolution  is:  1.13m 


3.21.2  Histogram 

Histogram  for  'Lwrun22-23'  Histogram  for  'Mwrun22-23' 


Figure  102  -  Gaussian  Fit  -  Histogram  of  Runs  22-23 


Histogram  Information  for  'Lwrun22-23' 

Normal  Model  Parameters: 

Mean  is:  1745.80 
Variance  is:  700.00 
RMS  error  is:  3.68% 

Chi2  P-value  is:  99.99% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Yes 
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Histogram  Information  for  'Mwrun22-23' 

Normal  Model  Parameters: 

Mean  is:  636.13 
Variance  is:  195.00 
RMS  error  is:  11.87% 

Chi2  P- value  is:  0.00% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  No 

3.21.3  Auto-covariance  -  Horizontal 

Auto-covariance  -  Horizontal  Slice  for  ’Lwnjn22-23'  Auto-covariance  -  Horizontal  Slice  for  'Mwrun22-23' 


Figure  103  -  Power  Law  Fit  -  Horizontal  Slice  of  Runs  22-23 

Auto-covariance  -  Horizontal  for  'Lwrun22-23' 

Mean  Standard  deviation  between  images  is:  36.44 
Power  Law  Model  Parameters: 
a_est  is:  941.40 
b_est  is:  0.00 
w_est  is:  0.12 
RMS  error  is:  4.31 
Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 

Auto-covariance  -  Horizontal  for  'Mwrun22-23' 

Mean  Standard  deviation  between  images  is:  21.43 
Power  Law  Model  Parameters: 
a_est  is:  360.60 
b_est  is:  0.00 
w_est  is:  0.09 
RMS  error  is:  6.33 
Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 
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3.21.4  Auto-covariance  -  Vertical 


Auto-covariance  -  Vertical  Slice  for  lLwrun22-23l 


Auto-covariance  -  Vertical  Slice  for  'Mwrun22-23' 


Figure  104  -  Power  Law  Fit  -  Vertical  Slice  of  Runs  22-2  3 


Auto-covariance  -  Vertical  for  ’Lwrun22-23’ 

Mean  Standard  deviation  between  images  is:  38.54 
Power  Law  Model  Parameters: 
a_est  is:  757.02 
b_est  is:  0.00 
w_est  is:  0.09 
RMS  error  is:  2.51% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Partially 

Auto-covariance  -  Vertical  for  ’Mwrun22-23' 

Mean  Standard  deviation  between  images  is:  22.22 
Power  Law  Model  Parameters: 
a_est  is:  357.42 
b_est  is:  0.00 
w_est  is:  0.09 
RMS  error  is:  3.95% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 
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3.21.5  Image  Synthesis: 


Sample  Ensemble  Image  for  'Lwajn22-23' 
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Figure  105  -  Original  Images  from  Runs  22-23 


Synthesized  Image  (Not  Transformed)  for  'Mwrun22-23' 
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Figure  106  -  Synthetic  Images  for  Runs  22-23 
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3.22  Run  24 

3.22.1  General  Information 

Flight  Date:  3-7-97DAY 
Altitude:  5K  ft 
Depression  Angle:  60  deg. 
Vegetation  Classes  in  Run:  18, 106 
Predominate  Vegetation  Class:  106 

Ensemble  Size  is:  350 
Number  of  Valid  Images  is:  174 
Vertical  Resolution  is:  0.65m 
Horizontal  Resolution  is:  0.56m 


3.22.2  Histogram 


Histogram  for  'Lwrun24' 


Greylevel 


Histogram  for  'Mwrun24‘ 


Figure  107  -  Gaussian  Fit  -  Histogram  of  Run  24 


Histogram  Information  for  'Lwrun24' 

Normal  Model  Parameters: 

Mean  is:  1753.09 
Variance  is:  490.00 
RMS  error  is:  4.05% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Partially 
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Histogram  Information  for  'Mwrun24' 

Normal  Model  Parameters: 

Mean  is:  696.84 
Variance  is:  126.00 
RMS  error  is:  4.24% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Partially 

3.22.3  Auto-covariance  -  Horizontal 


Auto-covariance  -  Horizontal  Slice  for  'Lwmn24'  Auto-covariance  -  Horizontal  Slice  for  'Mwrun24' 


Figure  108  -  Poxver  Law  Fit  -  Horizontal  Slice  of  Run  24 

Auto-covariance  -  Horizontal  for  'Lwrun24' 

Mean  Standard  deviation  between  images  is:  56.44 
Power  Law  Model  Parameters: 
a_est  is:  647.27 
b_estis:  1.31 
w_est  is:  0.11 
RMS  error  is:  4.23% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 

Auto-covariance  -  Horizontal  for  'Mwrun24' 

Mean  Standard  deviation  between  images  is:  15.36 
Power  Law  Model  Parameters: 
a_est  is:  148.29 
b_est  is:  0.00 
w_est  is:  0.11 
RMS  error  is:  3.96% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 
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3.22.4  Auto-covariance  -  Vertical 


Auto-covariance  -  Vertical  Slice  for  'Lwain24' 


180 

160 

140 

a 

O 

1 120 
(0 
s 

s100 

3 

< 

80 

60 

40 


Figure  109  -  Power  Lmv  Fit  -  Vertical  Slice  of  Run  24 


Auto-covariance  -  Vertical  for  'Lwrun24‘ 

Mean  Standard  deviation  between  images  is:  55.40 
Power  Law  Model  Parameters: 
a_est  is:  569.89 
b_est  is:  0.56 
w_est  is:  0.10 
RMS  error  is:  2.36% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Yes 


Auto-covariance  -  Vertical  for  'Mwrun24' 

Mean  Standard  deviation  between  images  is:  14.76 
Power  Law  Model  Parameters: 
a_est  is:  142.68 
b_est  is:  0.00 
w_est  is:  0.11 
RMS  error  is:  2.09% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Yes 


Auto-covariance  -  Vertical  Slice  for  'Mwrun24' 
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3.22.5  Image  Synthesis: 

Sample  Ensemble  Image  for  'Lwrun24' 
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Figure  110-  Original  Images  from  Run  24 

Synthesized  Image  (Transformed)  for  'Lwrun24' 
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Synthesized  Image  (Not  Transformed)  for'Lwrun24' 
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Figure  111  -  Synthetic  Images  for  Run  24 


Sample  Ensemble  Image  for  'Mwrun24' 
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Synthesized  Image  (Transformed)  for'Mwrun24' 
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Synthesized  Image  (Not  Transformed)  for  'Mwrun24‘ 


50  100  150  200  250 
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3.23  Runs  25-26 

3.23.1  General  Information 

Flight  Date:  3-7-97DAY 
Altitude:  5Kft 
Depression  Angle:  60  deg. 

Vegetation  Classes  in  Run:  106, 3, 47, 3, 105 
Predominate  Vegetation  Class:  3 

Ensemble  Size  is:  468 
Number  of  Valid  Images  is:  468 
Vertical  Resolution  is:  0.65m 
Horizontal  Resolution  is:  0.56m 


3.23.2  Histogram 


Histogram  for  'Lwrun25-26' 


Greylevel 


Histogram  for  'Mwrun25-26' 


Figure  112  -  Gaussian  Fit  -  Histogram  of  Runs  25-26 


Histogram  Information  for  'Lwrun25-26' 

Normal  Model  Parameters: 

Mean  is:  1732.32 
Variance  is:  460.00 
RMS  error  is:  4.98% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Yes 
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Histogram  Information  for  'Mwrun25-26' 

Normal  Model  Parameters: 

Mean  is:  690.32 
Variance  is:  118.00 
RMS  error  is:  5.28% 

Chi2  P-value  is:  100.00% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Yes 

3.23.3  Auto-covariance  -  Horizontal 

Auto-covariance  -  Horizontal  Slice  for  'Lwrun25-26'  Autocovariance  -  Horizontal  Slice  for  'Mwrun25-26' 


Figure  113  -  Power  Law  Fit  -  Horizontal  Slice  of  Runs  25-26 

Auto-covariance  -  Horizontal  for  'Lwrun25-26' 

Mean  Standard  deviation  between  images  is:  14.39 
Power  Law  Model  Parameters: 
a_est  is:  511.55 
b  est  is:  0.00 
w_est  is:  0.09 
RMS  error  is:  4.53% 

Model  Accepted  by  Threshold:  No 
Model  Accepted  Visually:  Partially 

Auto-covariance  -  Horizontal  for  'Mwrun25-26' 

Mean  Standard  deviation  between  images  is:  2.87 
Power  Law  Model  Parameters: 
a_est  is:  109.64 
b_est  is:  0.00 
w_est  is:  0.09 
RMS  error  is:  2.94% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Yes 
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3.23.4  Auto-covariance  -  Vertical 


Auto-covariance  -  Vertical  Slice  for  'Lwrun25-26'  Auto-covariance  -  Vertical  Slice  for  'Mwrun25-26' 


Figure  114-  Power  Lmv  Fit  -  Vertical  Slice  of  Runs  25-26 


Auto-covariance  -  Vertical  for  'Lwrun25-26' 

Mean  Standard  deviation  between  images  is:  15.04 
Power  Law  Model  Parameters: 
a_est  is:  503.56 
b_est  is:  0.00 
w_est  is:  0.09 
RMS  error  is:  1.61% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Yes 

Auto-covariance  -  Vertical  for  'Mwrun25-26' 

Mean  Standard  deviation  between  images  is:  2.99 
Power  Law  Model  Parameters: 
a_est  is:  112.62 
b_est  is:  0.00 
w_est  is:  0.10 
RMS  error  is:  2.25% 

Model  Accepted  by  Threshold:  Yes 
Model  Accepted  Visually:  Yes 
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3.23.5  Image  Synthesis: 


Figure  115  -  Original  Images  from  Runs  25-26 


Sample  Ensemble  Image  for  'Mwnjn25-26' 
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Synthesized  Image  (Transformed)  for  'Mwrun25-26' 
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Synthesized  Image  (Not  Transformed)  for 'Lwrun25-26' 
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Figure  116  -  Synthetic  Images  for  Runs  25-26 


Synthesized  Image  (Not  Transformed)  for  ’Mwrun25-2( 


100 


DSTO-TR-1456 


4.  Conclusions 


4.1  Histogram  Results 

As  the  analysis  of  this  data  shows,  the  majority  of  natural  terrain  background  is  non- 
Gaussian.  Figure  117  is  a  summary  of  the  results  showing  the  percentage  of  histograms 
which  passed  the  Gaussian  threshold,  which  were  a  partial  fit  and  which  didn't  pass 
the  hypothesis  test  (see  section  3.2  for  an  explanation  of  the  fit  criteria). 


Histogram  Summary 

60,  - - 


GaussianFit  Partial  Gaussian  Fit  No  Gaussian  Fit 


Figure  117  -  Histogram  Summary 

4.2  Auto-covariance  Results 

Not  all  of  the  runs  analysed  fitted  the  power-law  model.  To  demonstrate  the 
relationship  between  ensemble  size  and  standard  deviation.  Figure  118  and  Figure  119 
show  plots  of  the  standard  deviation  of  the  auto-covariance  slices  against  the  size  of 
each  ensemble.  They  indicate  that  as  the  ensemble  size  of  a  run  increases,  the  standard 
deviation  decreases  and  is  more  likely  to  fit  the  power-law  trend. 
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LWIR  Horizontal  Autocovariance  Slice 


LWIR  Vertical  Autocovariance  Slice 
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Figure  118  -  LWIR  Auto-covariance  Summary 

MWIR  Horizontal  Autocovariance  Slice 
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MWIR  Vertical  Autocovariance  Slice 
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Figure  119  -  MWIR  Auto-covariance  Summary 

4.3  Stochastic  Simulation  Results 


The  majority  of  the  transformed  synthetic  images  appear  to  visually  match  the  sample 
images.  However,  there  does  appear  to  be  some  differences  in  the  constrast  for  images 
which  are  not  transformed  but  pass  the  Gaussian  hypothesis  test.  This  difference  is  due 
to  distributions  having  pixel  values  that  are  either  really  small  or  really  large.  These 
cause  a  larger  dynamic  range  and  the  result  is  that  synthetic  images  have  an  incorrect 
constrast.  If  the  transformation  is  applied  to  these  images  before  synthesis,  they  are 
forced  to  be  truly  Gaussian  and  reasonable  simulated  images  result. 

The  more  structured  sparse  terrain  images  in  runs  10, 11, 13  and  15  also  do  not  appear 
to  match  the  real  images  closely.  Here  the  lightness  and  darkness  of  the  different 
regions  appears  correct  (a  property  of  the  histogram),  but  the  shapes  of  these  regions 
are  not  correctly  represented.  The  shapes  are  determined  by  higher  order  moments  of 
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the  random  field,  which  are  not  represented  by  the  histogram  and  auto-covariance 
estimates. 


4.4  General  Results  and  Further  Work 

As  the  1996  trial  suffered  from  high  atmospheric  absorption  and  low  temperature 
contrast  due  to  cloud  cover,  analysis  of  the  1998  trial  would  be  the  next  step  in  a 
complete  analysis  of  the  data  obtained.  Other  possible  avenues  of  work  arose  from  the 
following  observations: 

The  auto-covariance  power-law  model  may  not  be  the  best  fit  for  the  data  analysed  in 
this  report.  It  could  be  that  a  different,  more  complex  model  will  fit  the  data  and  help 
to  classify  the  correlations  in  natural  IR  imagery. 

The  non-Gaussian  nature  of  the  image  data  and  the  incorrectly  matched  shapes  for 
sparse  terrain  imagery  indicates  that  higher  order  statistics  are  important  for  a 
complete  statistical  analysis.  Chappie  et  al.  [3]  has  used  the  same  non-Gaussian 
simulation  model  and  compared  predicted  and  actual  third  order  autocovariance 
functions.  It  was  found  that  the  simulation  model  gave  a  reasonable  description  of  the 
third-order  statistics,  even  though  it  was  not  explicitly  incorporated  in  the  model. 

Alternative  approaches  to  the  stochastic  simulation  could  include  the  use  of  a  two- 
dimensional  Markov  field,  [15]. 
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